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Self-organized criticality in a discrete model for 
Smoluchowski’s equation with limited aggregations 

Mathieu Merle* Raoul Normand^ 

Abstract 

We introduce and study a discrete random model for Smoluchowski’s equation with 
limited aggregations. The latter is a model of coagulation introduced by Bertoin which 
may exhibit gelation. 

In our model, a large number of particles are initially given a prescribed number of 
arms. These arms are activated independently after exponential times and successively 
linked together. However, when the size of a cluster goes above a fixed threshold, it falls 
instantaneously into the gel, meaning that it no longer interacts with the other particles 
in solution. 

The concentrations in this model asymptotically obey Smoluchowski’s equation with 
limited aggregations. In this article, we study the discrete features of this model. We are 
able to argue that it remains closely related with a configuration model. We specifically 
obtain explicit expressions for the parameters of this configuration model, which show 
that it is subcritical before gelation, but remains critical afterwards. As a consequence, 
the asymptotic distribution of a typical cluster in solution is that of a subcritical Galton- 
Watson tree before gelation, while it is that of a critical Galton-Watson tree after gelation. 
Our model therefore exhibits self-organized criticality. 

Our study relies heavily on the study of the configuration model in the critical window. 
One of the main results of the paper is an extension of a result of Janson and Luczak on 
the size of the largest components in this critical window. 

Finally, as a consequence of our explicit expressions, we provide an explanation of ana¬ 
lytic formulas established by Normand and Zambotti regarding the limiting concentrations 
in Smoluchowski’s equation. 

Keywords: Gonfiguration model, gelation, hydrodynamic limit, limited aggregations, 
self-organized criticality, Smoluchowski’s coagulation equation 
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1 Introduction 

1.1 Smoluchowski’s equation 

The standard Smoluchowski’s equation was first introduced in 1916 in [30], to model 
coalescing particles. As time passes, clusters of particles are created. One assumes that 
each particle has unit mass, so a cluster has mass m, the number of particles it contains. 
As time passes, particles, and hence clusters, coalesce pairwise. The phenomenon is 
characterized by a symmetric kernel K{m,m'), modeling the “rate” at which two clusters 
of mass m and m' coalesce. When they do, a cluster of mass m + m' is formed. Then, 
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the evolution of the concentration ct{m) of clusters of mass m is given by Smoluchowski’s 
equation, actually an infinite system of nonlinear ODEs, 


—ct{m) = - E K{rn,m)ct{m)ct{rn — m)— ^ K{m,m)ct{m)ct{m), (1-1) 

^ / 1 /\1 
m —1 m'>l 


for m > 1. 

The behavior of this equation depends heavily on the kernel n. Three kernels are of 
particular interests, namely the constant, additive K,{m,m!) = m + m' and multiplicative 
K(m, m') = mm' kernels. For the first two kernels, the system is quite easy to solve, 
for instance by considering a PDE involving the generating function of {ct{m)). For the 
multiplicative kernel, three facts have been known for a long time. First, there is a unique 
solution on [0,r), where T = (X)m>i and it is mass-conservative, in that 

remains constant [22]. Secondly, if there is actually a solution on a larger 
interval, then the mass has to decrease after time T. This phenomenon is called gelation. 
Physically, it is interpreted as the appearance of infinite-size elusters accounting for a 
positive fraction of the total mass. Finally, Kokholm [20] has shown, for monodisperse 
initial conditions cfim) = that there exists a unique solution for all t >0, given 

explicitly for m > 1 by 


ct{m) 


m ^B{t, m) t < 1, 
{mt)~^B{l^m) t > 1, 


( 1 . 2 ) 


where 

B{\m) = 

m 


The function B{\,-) is a probability distribution for A < 1, called the Borel distribution. 
It is the law of the total progeny of a Gallon-Watson (GW) process with reproduction 
law Poisson(A), which we denote PGW{X). 

This shows in particular that the total mass J2m>i fncfim) is conserved up to time 1, 
but afterwards decreases and is equal to 1/t. Moreover, the proportion of particles which 
are in clusters of mass m is exactly the probability for a PGW{t) to have size m, when 
t < 1. On the other hand, when t > 1, the proportion of particles in clusters of mass m 
matches the probability for a PGW{1) to have size m. Since a PGW{\) is subcritical for 
A < 1, and critical for A = 1, this already suggests that the model exhibits a form of self- 
organized eritieality (SOG). Note also that the tail Z)m>A: ”^ct(m) decays exponentially 
with k before gelation, but as afterwards. 

The proof of existence and uniqueness to Smoluchowski’s equation with a multiplicative 
kernel and general initial eonditions has only recently been given in [26] , and independently 
in [27]. In the latter, it is also shown (for initial conditions with bounded support) that 
SOG occurs, in the sense described above that the tails of the total mass decay as 
after gelation. 

SOG is well-known in the Physics literature, but is seldom proved mathematically. 
Some models are however known rigorously to exhibit SOG, such as Abelian sandpiles 
[18], invasion percolation on regular trees [3], and Gurie-Weiss models [11]. 

Though the underlying idea in Smoluchowski’s equation is essentially discrete, all that 
we just said applies to the continuous equation, which is the natural mean-field model to 
consider. To obtain a discrete picture, we introduced in [23] a discrete model for Smolu¬ 
chowski’s equation, described as follows. Start with N particles, and give independent 
exponential clocks with parameter 1/N to each pair of particles. When a clock rings, form 
the corresponding link between the two particles, as long as they both belong to clusters 
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with sizes below a{N), where (a(iV)) is some threshold. This threshold should be thought 
of as when particles become too large and “fall into the gel”, that is, become inactive. This 
model is directly inspired by the work of Fournier and Laurengot [14], who first proved 
that, when the threshold satisfies 

1 < a{N) <C N, (1.3) 

then the concentration of particles of size m converges to the solution (1.2) to Smolu- 
chowski’s equation (1.1). To be precise, they only consider the concentrations as a jump 
Markov process, and prove convergence up to a subsequence, since their paper is an¬ 
terior to [26] and [27]. Their result is actually not limited to the multiplicative kernel 
K{m, m') = mm!^ but it is the only one which both exhibits gelation and is described by a 
truly discrete model as above. Rath in [27] studies a slightly different yet closely related 
model, and obtains the same convergence result. 

The threshold {a{N)) above which clusters become inactive is fundamental. Infor¬ 
mally, Smoluchowski’s equation (1.1) shows no interaction with clusters of infinite mass, 
so that it should also be the case in a relevant discrete model. In particular, there is also 
a convergence result for the more simple discrete model with no threshold {a{N)), i.e. 
when any cluster can interact with any other cluster, see [14] . The convergence is then to 
Flory’s equation, which is (1.1), but with an extra term taking into account interactions 
with clusters of infinite mass. As it stems from a more simple discrete model, this equation 
is also easier to study than (1.1), see e.g. [26] and references therein. 

Let us finally describe more precisely the results of [23]. The main interest of that 
model, compared to that of Fournier and Laurengot, is that it retains much more infor¬ 
mation, in particular the graph structure, if we think of particles as vertices and coagula¬ 
tions as edges being created. Therefore, it allows to study the shape of a typical cluster. 
Under stronger assumptions on the threshold, [23] shows that the shape a typical cluster 
converges to that of a PGW{t) tree for t <1, and PGW{1) for t > 1. In particular, this 
tree is subcritical before time 1, and critical afterwards, which is a far more striking way 
to explain SOC. As an aside, this also provides a better explanation for Formula (1.2), 
and a way to recover it. Indeed, the concentration of clusters of mass m is exactly the 
probability that a particle chosen uniformly at random is in an cluster of size m, times 
1/m, because of size-bias. From the result just mentioned, this probability converges to 
the probability that a PGW{t A 1) has size m, that is B(t A l,m), in agreement with 
Formula (1.2). See Sections 2.2.2 and 9.2 for similar reasonings and more details. 

1.2 Smoluchowski’s equation with limited aggregations 

Physically, one may want to think of a coagulation involving two given particles as the 
creation of a covalent bond, so that it makes sense to assume that each particle can a 
priori only form a certain finite number of such bonds. Alternatively, there could be 
several types of particles, some being able to create more bonds than others. Therefore, it 
seems reasonable to initially give each particle a certain number of arms, representing the 
number of covalent bonds it may create. Two arms are used to create a bond, or link, and 
cannot be reused again. Then, a cluster is characterized by its mass m and its number 
of free arms a, i.e. the total number of arms in that cluster which are not yet bound. 
This model with limited aggregations has first been considered by Bertoin in [5] , and gives 
the evolution of the concentration ct(a, m) of clusters with a free arms and mass m. It is 
given by an equation resembling Smoluchowski’s equation (1.1), which we now describe. 

Let 5 = N X N*, and for p = {a,m),p' = {a',m') e 5, define p ■ p' = aa' and 
pop' = (a -|- a' — 2, m -|- m'). We write p' <p a' < a+ 1 and m' <m — l. For p' < p, we 
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define p\y = (a + 2 —a',m —m'). Then Smoluchowski’s equation with limited aggregations 
is given for p G S' by 


dt 


ct{p) 


^ ^ p' • {p\p')ct{p')ct{p\p') - p-p'ct{p)ct{p'). 

p'<p p'&S 


(1.4) 


The first term in the RHS accounts for the appearance of (a, m)-clusters, by coagulation 
of (o', m')-clusters with (a + 2 — a',m — m')-clusters. The second term accounts for the 
disappearance of (a, m)-clusters by coagulation with any other cluster. 

In [5], Bertoin considers monodisperse initial conditions co{a,m) = for a 

measure^ p! on N with unit mean and a second moment, and such that p'{d) + p'{2) ^ 1. 
Denote 7 ' = X)a>i “ 2)/i'(a). Bertoin then shows that (1.4) has a unique solution up 
to time Tggj, interpreted as the gelation time, with 


0 ^ / + 00 , 

U/v. 


if 7 < 0 
if 7 > 0 . 


(1.5) 


More recently, Normand and Zambotti have shown in [26], using analytic techniques, that 
for general initial conditions (and in particular even when 7 ' > 0), Equation (1.4) has a 
unique solution on the whole of M"*". They also provide explicit solutions; we will come 
back to this shortly. 

One very interesting formula of [5] is the following. Denote v{m) = (m + l)p'{m + 1), 
which is a probability measure on N. Then, provided 7' < 0, there are limits to the 
concentrations, given for m >2 and a > 0 by 

Coo(a,m) = ^ Y ^im-2), (1.6) 

where is the m-th convolution product of u. One part of this formula is clear: at the 
end, all the arms have been used up. To provide an explanation for Coo(0,m), Bertoin 
and Sidoravicius [ 6 ] study a microscopic model for these limiting concentrations. Without 
dwelling on technical details, they consider N particles, with an empirical distribution of 
arms converging (weakly and in mean) to p'/p'(N). They then consider a uniform 
pairing of these arms. Then, in the case when 7 ^ < 0, they show that, when an arm is 
selected uniformly at random, the distribution of the cluster to which it belongs converges 
to that of GW tree started from 2 individuals, and with reproduction law n. This allows 
to give a probabilistic interpretation of Formula (1.6), see [ 6 ]. We will provide another 
interpretation of (1.6) in Section 2.2.2, based on the configuration model and which is 
more in the flavor of this work. 

Let us now explain what happens when gelation does occurs, i.e. in the case when 
7 ' > 0. Write for the generating function of n, and assume z^(0) > 0. Then it is easy 
to see that the equation 

cGUc) = G,{c) (1.7) 

has a unique solution c G (0,1), and that /3 = 1/G(^(c) > 1. In [26], it is then proved that, 
for all m > 2 and a > 0 , 

Coo(a, m) = l{a=o} , ^ - 2). (1.8) 

‘ m[m — Ij 

^The prime is to differentiate Bertoin’s convention and ours, where he normalizes p to have unit moment, 
whereas we normalize it to be a probability, and we will then drop the prime. 


4 







This is obviously extremely similar to (1.6), up to the factor which seems to favor 

larger clusters. A question left open in [26] is to give a probabilistic interpretation for this 
formula. We propose to answer this question, and more importantly, study in details a 
microscopic model for (1.4), which, similarly to the model of [23], exhibits SOC. 

A final result in [6] is a dynamic model for (1.4), where each pair of free arms is linked 
at rate 1/N. The authors then show that the concentrations of clusters with mass m 
and a free arms at time t converges to the solution ct{a,m) of (1.4). Of course, its final 
state corresponds to a uniform pairing of the arms, in accordance with the reasonings 
above. This is the natural discrete model underlying (1.4), but, as should be clear from 
the previous section, it is not quite the right model when gelation occurs. 

1.3 Model 

Considering that we already explained the models of [6] and [23], the model that we are 
now going to present should not come as a surprise. So let us start initially with N 
particles, each being equipped with a certain (deterministic) number of arms. We denote 
the degree distribution, that is, there are N particles with r arms for every 
r G N. We also fix a threshold a{N). As time goes by, arms will bind to create clusters. 

We shall call a cluster small or in solution if it has size less than a{N), and large^ or 
in the gel if it has size a{N) or greater. Similarly, we will say that a particle or an arm 
is in solution if it belongs to a small cluster, and otherwise that it is in the gel. We also 
say that an arm is free if it is not linked to another arm, so that they are all free at time 
0 . 

Now, each arm is given independently an exponential clock with parameter one £(1). 

When a clock rings, the corresponding arm is made active. If, at this instant, there is 
another active and free arm in solution, the two arms are linked together to form a bond 
between the corresponding particles (and are then no more free). Otherwise, nothing else 
happens. 

Note that there can only be zero or one active arm in solution. Another way to see this 
model is to wait for an arm to activate, then for the next one, and bind them. Continue 
with the third and fourth activated arms, and so on, until gelation occurs, that is, a large 
cluster is created. This cluster then precipitates, or “falls into the gel”, which is just a 
metaphor to mean that it becomes inactive, and it cannot bind to any more particle. The 
process continues similarly on the particles in solution, until another cluster falls into the 
gel, and so on. 

The reader might notice two differences with the dynamic version of the model of 
Bertoin and Sidoravicius, which we explained above. First, of course, there is the appear¬ 
ance of the threshold ct{N), which is fundamental, as for (1.1). Without it we would get 
convergence to the corresponding Flory equation discussed in [26]. On the other hand, 
we put clocks on arms, not on pair of arms. It is clear that the only difference it creates 
is a change of time, but the underlying Markov chain is the same, since, in both cases, 
a coagulation occurs by picking a pair of free arms in solution uniformly at random. We 
choose this model because the microscopic evolution is more clear. Indeed, to know what 
happens when a clock rings, we just need to check the state of one arm. Otherwise, 
we need to check the state of both arms in the link, which leads to some cumbersome 
complications. As expected, and as we shall see in Section 3, the convergence of the con¬ 
centrations will be to a modified Smoluchwoski equation with limited aggregations, which 

^We use large in contrast to giant, which is the term usually used in the Random Graphs literature to describe 
clusters of size of order N. 
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is just a time-change of (1.4). 


1.4 Configuration model 

Consider the easier version of our model, where we also set i.i.d. £(1) clocks on the arms, 
but there is no threshold a{N). Hence, an arm is activated when the corresponding clock 
rings, and at this time, if there exists another free active arm, we bind these two arms 
together. In other words, we bind the first and second activated arms, the third and fourth 
activated ones, and so on, up to time t. As we mentioned, this is not the right model for 
(1.4), but it is nonetheless relevant, and it is also easier to study. 

At the end of this process t = -|-oo, we have thus bound all the arms (except maybe 
one), to obtain a random (multi-)graph. Notice that it would be exactly the same, in law, 
to successively choose a uniform pair of free arms until they are all picked, or to choose 
a pairing of the arms, uniform among all the possible pairings. This is the well-known 
configuration model (CM) with degree distribution see e.g. [9, 29]. 

In this modified model, if we stop at some time t, then each arm is activated inde¬ 
pendently with probability 1 — e“*, but, conditionally on the activated arms, the order 
of activation is uniform. We thus still deal with a CM, but whose degree distribution 
is random. Note however that, by the law of large numbers, this degree distribution is 
essentially constant. The process described above will thus be called a dynamical CM 
(DCM). 

There is an obvious coupling between our model and a DCM, by using the same 
exponential clocks, and these two processes are exactly equal up to gelation, that is, the 
time when a large cluster appears. Therefore, it seems reasonable that information on 
(dynamical) CM would give information on our model. For instance, if we know the size 
of the largest component in a DCM, then we know the gelation time in our model. This 
is described in more details in Section 2.3. 

Explaining our results takes same time and preparation, and we shall delay this to 
Section 2. We first conclude this introduction with a short review of the literature, and 
the description of the plan of our paper. 

1.5 Related works 

1.5.1 Frozen percolation on trees 

Aldous considers in [2] a model of frozen percolation on a binary tree. Edges are given 
independent uniform on [0,1] clocks, and when a clock rings, the edge appears if and only 
if both vertices belong to finite clusters. A rewording of the latter condition is that any 
cluster which becomes of infinite size is instantaneously frozen and can no longer interact 
with other clusters. 

Gelation happens when the first infinite cluster is formed, at time 1/2. After proving 
existence of such process past 1/2, Aldous computes in particular the rate at which the 
cluster containing a given edge or a given vertex becomes frozen. He further shows that 
at any time past gelation, conditionally given that an edge e is present and its cluster 
is finite, clusters of the edge e is a critical Galton-Watson tree with two ancestors. On 
the other hand, the distribution of infinite clusters is always that of an incipient infinite 
cluster (in this case a Galton-Watson with Bin(2,1/2) offspring distribution conditioned 
to be infinite). 
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1.5.2 Forest-fire models 

Rath and Toth in [28] prove SOC for a model of forest fires on Erdds-Renyi random 
graphs. In their model, edges between each pair of the N particles appear (or re-appear) 
at rate 1/A^, while lightning strikes particles at a fixed rate a{N)/N, with {a{N)) satis¬ 
fying Assumption (1.3) for the most interesting behavior. The effect of lightning is that 
edges in the corresponding connected components simply vanish (but vertices remain). 
The characterization of SOC in [28] is through the fact that past gelation, the limiting 
concentrations {ct{m), m > 1) satisfy J2m>k mctim) C{t)k with C'(t) > 0 for t > 1; 
for t < 1, the decay is at an exponential rate. 

Related to this model. Rath considers in [27] a model of frozen percolation on the 
Erdos-Renyi random graph. Edges between each pair of the N particles appear at rate 
1/N, and vertices are hit at rate a{N)/N. The effect of a vertex being hit is that its 
connected component becomes instantaneously frozen. Rath obtains convergence of the 
concentration of clusters of size m towards the solution to Smoluchowski’s equation with 
a multiplicative kernel. He further proves SOC for this model, in the same sense as [28] 
above, and obtains limit theorems for the remaining mass and the distribution of the size 
of the component of a typical vertex. 

Although the models of [23] and [27] present obvious similarities, they are also quite 
different: for instance the concentrations in the model of [23] converge to the solution 
of Smoluchowski’s equation with a multiplicative kernel, but the ones in [28] do not. 
Intuitively, this is due to the presence of additional dust created by the burning of clusters, 
which then also affects the growth of larger clusters. The characterization of SOC in these 
two papers is also quite different, since [27] obtains SOC on the continuous model, while 
our goal in [23], as in the present paper, is to observe SOC on a microscopic level. 

1.6 Organization of the paper 

Following this. Section 2 is a precise outline of the paper, showing rigorous statements 
of our results and the main steps and ideas of the reasoning. We also introduce our 
assumptions and the main notation. For convenience, an index of notation is added in 
Appendix A. Section 2 contains in particular an explanation of the link between our model 
and the configuration model, and gives the main idea of how the latter helps study the 
former. 

In Section 3, we state and prove the main results concerning Equation (2.8). This 
section is essentially independent of the rest of the paper, and is merely meant to justify 
that we indeed get a model for Smoluchowski’s equation with limited aggregations. In 
particular, the details can be skipped in a first reading. 

We get back to our model in Section 4. This is a part containing three types of 
preliminary results. First, that, for any time, there is a still a positive proportion of 
particles remaining in solution, what allows to use asymptotic results. Secondly, some 
fairly easy technical results which allow, thanks to our assumptions, to strengthen most 
convergence results that we get. Finally, we prove some combinatorial results describing 
the structure of our system at any time in terms of the configuration model, what allows 
to justify the comparisons between the two models. 

The longest part of this article is by far Section 5, where we study the CM in the 
critical window, and prove an improved version of the result of Janson and Luczak. That 
section is focused on the CM and is entirely self-contained. It has therefore an interest of 
its own for readers who are more interested in random graphs. We also prove an extension 
of this result for the percolation on a CM. 
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These results are applied to dynamical CM in Section 6, specifically to study when, 
in a DCM, we observe the appearance of a cluster of size a{N). We also determine 
the properties of this cluster. This is a quite straightforward application of the previous 
section, but the results are slightly tedious to write, and require a lot of notation. 

To clarify the link between our model and DCM over several gelation events, we 
introduce in Section 7 an alternative model. It is built from independent DCM, more 
precisely from parts of DCM pasted together at the right times. We construct in this way 
a coupling of our model with another easier model, and we prove, using the results of the 
preceding section, that they coincide with high probability, at least on compact intervals. 

This completes the most technical part of the work. What remains is to use this 
comparison, and the results on the DCM, to study our process. This is done in Section 8. 
We first prove tightness of our processes, and then show that subsequential limits have to 
solve a PDE. Despite the unusual nature of this PDE, it turns out to be exactly solvable. 
This allows to conclude in a classical fashion. 

Einally, Section 9 is a conclusion, where we first explain how to recover (1.8), and 
solve our initial problem. We also present some examples of initial distribution for which 
computations can be carried out explicitly. We conclude by discussing our assumptions 
and some open problems. 

One last note to the reader: [23] and the present paper obviously use similar ideas, 
but the presence of arms creates many complications, both in the reasoning and in the 
techniques (as the length of the papers can attest). We thus think that [23] is a good 
preliminary reading. 


2 Outline and main results 

The goal of this section is to introduce some notation, describe precisely our results, and 
finally get a flavor of the reasoning route that we will follow, by studying what happens 
when gelation occurs. 

2.1 Notation and assumption 

Eor convenience, we summarized all the notation that we use in Appendix A. To begin 
with, for a measure /r on N = {0,1,2,... } and i e N, we let 

G^{x) = xe[0,1], 

k>0 

for its generating function, and 

mf = G«(l) = J2Hk-l)...{k-i + l)f,{k) 
k>0 

for the i-th factorial moment of /i. Weak convergence of measures is denoted by —?■. 

Recall that in our model, we start with N particles indexed by [N] = {1,2,...,A^}, and 
with degree distribution . Clusters of size < a{N) are small (or in solution), those 
of size > a{N) are large (or in the gel). We set i.i.d. exponential clocks with parameter 
one on each arm, and bind successively the first and second activated in solution, the 
third and fourth in solution, and so on. We let G{t) describe the graph created by the 
particles and their bonds at time t. In fact, we will always consider graphs (and trees) up 
to graph isomorphism, but this is a minor technical issue that we will not mention again. 
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For the sake of definiteness, we assume all our processes to be cadlag. When we talk 
of the degree of a particle, we always mean its total number of arms (activated or not), 
which is therefore a constant in time for a given particle. 

We will only be interested in particles in solution, so we let Gs{t) be the restriction of 
G{t) to the particles in solution at time t. Note that there are no bounds between Gs{t) 
and its complement. We also define 

• Nt, the number of particles in solution; 

• nt = Nt/N, the concentration of particles in solution; 

• Pt{k,r), the number of particles in solution with k activated arms and degree r; 

• pt{k,r) = Pt{k,r)/N, the concentration of particles that are in solution and with k 
activated arms and degree r; 

• /Ut(r) = J2k=o Pt{k,r)/Nt, the distribution of the degree of the particles in solution; 

• TTt{k) = ^t{k,r)/Nt, the distribution of the number of activated arms of the 

particles in solution. 

With a slight abuse of notation, we have po = /r. Note that all these quantities are 
encoded in (pt), and that pt and nt are probabilities. When a large cluster appears, we 
witness a gelation event. The first one will hold a special role and will just be called 
gelation. We will denote 

• the time of the successive gelation events; 

• "^gei = Ti, the gelation time. 

When necessary, we will add an exponent (iV) to these quantities to mean that we talk 
of the model on N particles. Finally, we say that a sequence of events occurs with 

high probability (w.h.p.) if P{E]\f) —)■ 1 as N — +oo. 

We will need two types of assumptions. One concerns the degree distribution 
the other the threshold (a(N)). 

Assumption 2.1. There exists a probability p sueh that p^^'^ —>■ p and > 0. Moreover, 
there are eonstants rj,M ^ (0, +oo) sueh that 

sup Y^p^^\r)r^+^ < M. 

A^>lr>0 

„(iV) 

Notice that the second part of the assumption also ensures that m/ —)■ m/ for 
i = 1,2, 3,4. It turns out that this assumption will suffice for results on the CM, but that 
we will need a bit more so as to study dynamieal CM, and our own model. 

Assumption 2.2. There exists a probability p ^ do sueh that p^^'l —)• p and rrig > 0. 
Moreover, there are eonstants rj,M ^ (0, +oo) sueh that 

sup y^p^^\r)r^+^ < M. 

The only difference is of course the 5+r/-th moment assumption. Now, as we mentioned, 
the result of convergence to Smoluchowski’s equation only requires a simple assumption 
on (a{N)), namely (1.3). As in [23], we will need more on (a{N)) to be able to use results 
on the configuration model, and actually even more than the “minimal assumption” 

< a{N) <C N (2.1) 

so as to be able to use CM results repeatedly. 
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Assumption 2.3. As N ^ +oo 


<C a{N) < N. (2.2) 

The careful reader will notice in the proofs that there is a balance between the as¬ 
sumptions on and on {a{N)), in that stronger moment conditions allow to take 

smaller a{N). We did not try to optimize this. 


2.2 Configuration model 

2.2.1 Phase transition 

Consider N vertices, with N7r^^\i) of them of degree i (that is, endowed with i arms, 
or half-edges). The probability is called the degree distribution. Performing a 
uniform pairing of all the arms (with eventually one half-edge left unpaired, if their total 
number is odd) gives a random graph, called the configuration model 
Let us denote by Ci(A, and Q 2 {N, vr^'^^) the largest and second largest components 
of the graph CM{N, (chosen uniformly at random if there are several choices). We 

write and 62 ^^ if there are no ambiguities. For a given component C, we let Vk{G) 
be the number of vertices of degree k in S, and u(S) = Z)fc>o'^fc(C) its size. 

The configuration model has been extensively studied since its introduction by Bollobas 
in [ 8 ] . There are in fact several variants for its definition, see for instance Chapters 7 and 
10 of [29] for a survey. In this paper we take the simplest definition as a uniformly chosen 
pairing of arms, thus allowing the possible presence of self-loops or multiple edges, and 
leading to an obvious (multi-)graph structure. We also allow one arm to be left unpaired 
when there is an odd total number of arms, without affecting the graph structure. 

Denote and let us assume that converges weakly to some proba¬ 
bility TT, with —)• mj and —)• mj. An important parameter for the CM is 


7 


TT 


777,2 ~ • 


(2.3) 


An abrupt transition between the subcritical regime < 1, where most vertices are 
in finite components, and a supercritical regime 7 ^ > 1 , where a positive proportion of 
vertices belongs to a single giant component, was shown by Molloy and Reed in 1995 [25], 
up to some technical conditions, which were removed in [16]. We will be more precise 
shortly, but let us first start by explaining what happens at a microscopic level. This 
should enlighten the appearance of this parameter 7 , as well as the following results. 

For a measure vr 7 ^ do on N, we will denote vr the size-biased, shifted by one 
probability defined by 


7t{k) 


{k + l)-n{k + 1 ) 
Ei>i*7r(i) 


fe e N. 


(2.4) 


To see why this is a relevant notion, look at a CM(N, vr) for large N. One way to build 
it is to select a vertex uniformly at random, take one of its half-edges e (chosen in any 
way), then choose another half-edge ^ uniformly at random, and create the link between 
e and e^ Then, select another free half-edge in this cluster (in any way), choose another 
free half-edge uniformly at random, bind them, and so on. Eventually, we may run out 
of half-edges in the cluster, and we may then start with another vertex. In this way, we 
then build a CM component-wise. This construction is actually at the heart of the proof 
of [16], and will be used in Section 5. 

Now, let us see what happens in this algorithm. First, we pick a vertex v uniformly at 
random. Its number of half-edges has law ~ tt. Then, we pick an half-edge uniformly 
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at random, belonging to a vertex v'. Forgetting about the half-edge we already chose, it is 
easy to see that the number of half-edges of v' is distributed according to ~ vr. The 
difference clearly lies in the fact that we choose a half-edge, rather than choosing a vertex. 
Repeating the result, we see that, as long as the cluster is not too big, we essentially 
always pick new vertices, and choose them according to vr. Therefore, the cluster we will 
see will be a tree, whose root has a number of children distributed according to tt, and 
whose other vertices have a number of children distributed according to tt. In branching 
processes terminology, this is a delayed Galton-Watson tree, that we denote 

So what does it have to do with 7 ? Recall [4] that a Galton-Watson tree GW,r with 
reproduction law vr has extinction probability one if and only if rnl < 1. Therefore, it 
is easy to see that a GWjr,* has extinction probability one if and only if rrii < 1. But 
one readily computes and therefore there is a.s. extinction if and only if 

m^lrnl < 1, that is, 7 ^ < 0. Though this is informal, this explains the appearance of this 
quantity For instance, a similar reasoning could be done for the Erdds-Renyi graph. 

In the next section, we will mention some results showing an even more striking relation 
with Galton-Watson trees, but let us first state the phase transition result. 

Theorem 2.4 ([15, 16, 25]). Assume that — )■ vr with 7 r(l) > 0 and 

sup ^ TT^^\k)k‘^ < -|-oo. 
k>0 

Then 

1- if 1 '" < 0 , then v{q[^'^)/N —)• 0 in probability; 

2. if> 0, then v{Qi^'^)/N converges in probability to the (positive) survival proba¬ 
bility of a GWTT^ir-tree. 

If moreover 

sup "S) 'K^^\k)k‘^~^^ < -|-oo 
fe >0 

and 7 ^ < 0, then there is a constant K such that w.h.p. 

In other words, a giant component exist with positive probability if and only if 7 ^ > 0. 
If it does, its size is that of the survival probability of a GW 7 r, 7 T-tree, similarly to the 
Erdos-Renyi case [17, 29]. This is the classical result of [25], with the assumptions of [16]. 
These two papers do not make precise the size of the components when 7 ^ < 0. This 
is however done in [15], under some growth assumptions. Under the assumption done, 
corresponding to Assumption 2.1, it amounts to the statement above. 

In parallel to the Galton-Watson terminology, we will therefore call the cases 7 ^ < 0, 
7 ^^ = 0 and 7 ^^ > 0 respectively subcritical, critical and supercritical. Notice in particular 
that Tggj, defined in (1.5), is finite if and only if we are in the supercritical case. We will 
explain this in section 2.3. 

2.2.2 Local convergence 

We explained above what a typical component in a CM{N, vr^^^) should look like, namely 
a This intuition can actually be made precise through the notion of local con¬ 

vergence'. essentially, a rooted graph converges locally to another rooted graph if they 
look the same at any fixed distance around the root. Eor probabilistic results, the natural 
extension is the concept of local weak convergence popularized by Aldous, see e.g [I, 10]. 
The following result is exactly Theorem 3.15 of [10], see also Proposition 2.5 of [12]. 
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Proposition 2.5. [10, 12] Assume that —)■ vr with —?■ mj. Consider a vertex 

uniformly at random, and the eomponent rooted at this vertex. Then, in the sense 
of loeal weak eonvergenee 

as N ^ +00. 

This makes the intuition from the previous section rigorous: a typical component 
of a CM asymptotically looks like a delayed GW tree. In particular, in the subcritical 
case, these trees are finite a.s., and the distribution (as a graph) of a typical component 
converges to that of a tree. 

This result of weak local convergence allows to give a more direct explanation to this 
formula than the one given in [6] that we mentioned in the introduction. Indeed, c^^(0, m) 
is the concentration of clusters of size m (let us forget about the possible unique free arm 
remaining) in a CM{N There is a size bias when one picks a vertex uniformly at 
random, and thus 

where Q is the probability corresponding to the choice of the root vertex of (so the 

RHS above is still random). But under the assumptions of Proposition 2.5, —)■ GWjr^*, 

and therefore 

Q (u(e(^^) =m) (u(GW^,*) = m) = —3_7r*™(m - 2), 

what explains Formula (1.6). The last part of the above computation is a simple applica¬ 
tion of the celebrated formula of Dwass [13] . 

We will use the same type of reasoning to explain Formula (1.8) in Section 9.2. How¬ 
ever, in this case, we do not know what a typical cluster look like. This is what we will 
investigate in this article, as explained in the rest of this outline. 


2.2.3 The critical window 


The reader will have noticed that Theorem 2.4 does not provide much information about 
the case = 0. Just as in the case of the Erdds-Renyi graph, what happens in this 
eritieal window is more subtle. So let us first state the result of [16] describing precisely 
what happens then. We let y^v = 

Theorem 2.6 ([16]). Suppose that Assumption 2.1 holds for vr. Assume moreover that, 
for jN = , 

<C yAT < 1. 

Then, for every J > 0, the following hold w.h.p. 

1. The size of the largest eomponent verifies 


V 




< JiVyAT. 


2. For every k >0, 



- 2 


k'K{k) 


nin 


iVyAT 


< dNjisf. 


12 








3. The size of the seeond largest eomponent verifies 

V (ef)) < 5N^n. 

This result describes precisely what happens in the critical window. First, Points 1 
and 3 show that there is a unique largest component, which has size about 
Notice that this quantity is of an order strictly between and N^ which is the reason 
for Assumption 2.1. The second point means that the vertices in that large component are 
chosen by a size-bias of their degree. This amounts to saying that, to find the vertices in 
that component, one chooses about /rri^ arms uniformly at random, and picks 

the corresponding vertices. 


2.3 First gelation time 


Let us explain how to use these results to study our model. We assume that the initial 
distribution of degrees verifies Assumption 2.1 and that {a{N)) verifies (2.1). Recall 

that G{t) is the graph induced by the particles at time t, along with their links. As we 
already mentioned in the introduction, there is another model coupled with ours that we 
can consider: link the first and second activated links, the third and fourth, and so on, 
until time t, to get a graph process (S(t))- In this other process, that we call a dynamic 
configuration model (DCM), there is no gel and solution, and any cluster can interact 
with any other cluster. 

Recall that Tgei is the gelation time, that is, when we first see a large component. It is 
clear that 

G{t) = 9{t), t< Tgel. 

In particular, Tgei is exactly the time when a large component appears in the DCM 9, and 
the size and properties of that component are the same for both models. Therefore, up 
to and ineluding Tgei, it suffices to study 9- 

Now, for some fixed time t, 9(t) is a CM on N vertices, with a random degree distribu¬ 
tion p\^\ where corresponds to picking an arm or not with probability pt = 1 — e“*. 
By law of large numbers, it is easy to check that p\^'^ —>■ pt a.s. with 

Gpfix) = Gp ((1 - e“*)x -h e“^) , 

and even more, than Assumption 2.1 holds a.s. for {p\^'^). The generating function allows 
to compute 

ml* = p\ mf, 7^* = Pt {pt mif - mf). 

Clearly, is an increasing function of t that vanishes at Tgei whenever it is finite, where 


Tgei — 



if > mf 
otherwise. 


(2.5) 


Notice that Tgei is finite if and only if T^gj, defined in (1.5) is finite. We can then see 
different alternatives. 


The subcritical case The subcritical case corresponds to 7^ = — mf < 0 and 

is simple. Indeed, when all the arms are activated, the graph 9(oo) is a GM{N, 
which is subcritical. In particular. Theorem 2.4 implies that its largest component has 
size 0 {N^^G+v)'j ot{N), since we assume (2.1). This means that w.h.p. there is no 
gelation event, and the final configuration of clusters is that of this configuration model. 
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The critical case The critical case corresponds to 7 ^ = 0, and is not precise enough 
to decide whether or not there will be at least one gelation event (for that, one would 
at least need to know the leading order of 77 V as N ^ 00 , according to Theorem 2.6). 
Nonetheless, at any time t >0, 7 ^‘ < 0, and therefore the graph S{t) is an asymptotically 
subcritical CM. By the exact same reasoning as above, w.h.p. there is no gelation event 
up to time t, and the configuration of clusters at any time t > 0 is that of a configuration 
model with N particles and degree distribution pt. 

Notice that the critical and subcritical cases mean that 

• ^gel — TCX), 

• Tggj = + 00 , that is, there is no gelation in Smoluchowski’s equation (1.4); 

• there is w.h.p. no gelation (at a finite time) in our discrete model. 

In other words, —)• Tgei = + 00 . As expected, when Tgei < + 00 , that is, there is 
gelation in ( 2 . 8 ), we will observe gelation in our discrete model as well, as we now explain. 


The supercritical case The supercritical case corresponds to 7 ^ > 0. Then 7 ^* < 0 
for all t < Tgei, and by a reasoning similar as above, there is w.h.p. no gelation event 
at a time t < Tgei. On the other hand, having no gelation event by time t > Tgei means 
that a CM{N, p\^^) has no cluster of size greater than a{N) <C N. But then 7 ^* > 0, so 
by Theorem 2.4, this event has asymptotically vanishing probability. Hence, this simple 
argument tells us that 

iT ^ T^gel, (2.6) 

provided again that (2.1) holds. Clearly, this even just requires that a{N) <C 

N. 

Using Janson and Luczak’s result in the critical window. Theorem 2.6, we can in fact 
be even more precise in understanding typically when the first gelation event occurs, and 
what happens then. Let us indeed consider the time 


tN 


a{N) mg 
N 2m2{m2 — rrii) 


Since a{N)/N —>■ 0, it is easy to check that 


In 


(N) 

P) 


a{N) 

N 




Now, {a{N)) verifies (2.1), and thus the assumptions of Theorem 2.6 are verified. The 
largest component in a CM {N, ) has then size about 

711 ~ 

2iV—^7Ar ~ a{N). 

Therefore 

_^(7V) ^ ^3 

N 21712 {m^ — rrii) \ A" / 

Moreover, we see that the component that is created at this time T^^^ has size 
a{N) + o{a{N)), when we can a priori only tell that it is between a{N) and 
Theorem 2.6 even allows to be more precise, and to show that, in that first gel 
the particles are chosen by a size-bias of their number of activated arms. 


(2.7) 

exactly 

2a{N). 

cluster. 
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Of course, we are slightly cheating here; the result (2.6) holds rigorously, but (2.7) is a 
bit informal. To begin with, we did not make explicit what o(-) means. We also dealt quite 
liberally with the error terms. This could be hxed easily by taking times slightly before 
and after tn, as we will do in Section 6. A most important imprecision however, is that we 
used reasonings on the CM with a random degree distribution, without checking that the 
assumptions of Theorem 2.6 hold a.s. with respect to the law of the degree distribution. 
As will turn out, this will require slightly more to be made precise, to wit Assumption 
2.2. We refer to Section 6 for details. 


2.4 Results 

As should be clear from the earlier description, there are three main protagonists in our 
play: Smoluchowski’s equation with limited aggregations (1.4); our discrete model; and 
the conhguration model. We now describe our results concerning these objects. 


2.4.1 Modified Smoluchowski equation 

To begin with, we will be concerned with Smoluchowski’s equation with limited aggrega¬ 
tions (1.4), or rather the slight modihcation that we now introduce. For / : 5 —)■ M"*", we 
write {ct, f{a,m)) = Z)(a,m)eS/(“’™') equation is 


dt 


ct{p) 


1 

2{ct,a) 


^ p • p\p'ct{p)ct{p\p) - pct{p)ct{p) , 

\p'<p p'&S J 


( 2 . 8 ) 


for p € S. Two natural questions are then raised. First, is this equation well-posed? 
Secondly, what is its relevance to our model? The two following answers should not come 
as a surprise. 

Let us consider initial conditions cq, with 


nil = (co, a) G (0, -|-oo), m 2 = (cq, a{a - 1)) G [0, -|-oo] 


and let 


Tgel — 


- log (1 

-|-oo 


if "^2 > mi 
otherwise. 


(2.9) 


Clearly, for monodisperse initial conditions co{a,m) = /r(a)l|m=i} for a probability p on 
N, then this is the same as the time Tgel dehned in (2.5). Since, we will always consider 
monodisperse initial conditions, we keep the same notation. The hrst result is the weh- 
posedness of (2.8). The precise dehnition of a solution is given in Section 3. 


Theorem 2.7. Consider initial conditions {co{p),p G S) with (co,a) G (0,-|-oo), and de¬ 
fine Tgel as in (2.5). Then the modified Smoluchowski’s equation with limited aggregations 
(2.8) has a unique solution c on M^. It enjoys the following properties. 

1. The mass {ct,m) is continuous on M+, constant before time Tgel, and decreasing 
afterwards. In other words, there is gelation at time Tg^i. 

2. The total number of arms At = {ct,a) is continuous, decreasing and remains positive 
on M"''. 

3. The ODE s'{t) = A^(j),s(0) = 0, has a unique solution on M_|_. This solution is a 
C^-diffeomorphism ofMA and is the solution to (1.4) with initial condition cq. 
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4- If b is the solution to (1.4) with initial conditions co, then for Bt := {bt,a), the 
ODE s'{t) = 'S(O) = 0, has a unique solution on M"''. This solution is a 

-diffeomorphism o/M"*" and bg(i^ is the solution to (2.8) with initial conditions cq. 

5. The limiting concentrations Cooip) = q(j>) exist for all p G S and are the 

same as for (1.4), that is, given by (1.8). 

This shows in particular that Equations (1.4) and (2.8) are essentially the same, since 
they are related by a time-change. Now, as we mentioned. Equation (1.4) corresponds 
informally to a model where pairs of arms are bound at rate 1. Since in our model, we 
have to wait for two activations to bind a pair of arms, we should observe a time-change at 
the limit, which turns out to give Equation (2.8) instead. Indeed, if we denote c\^\a,m) 
for the concentration of clusters with a free arms and mass m in our model, then we have 
a result similar to that of [14], provided we also assume (1.3). As usual, ©(M"*",!!!) is the 
Skorokhod space of cadlag functions on M"*" with values in a Polish space E, see [7, 19] . 

Theorem 2.8. Assume that (1.3) holds and that converges weakly to a probability 
pL, with 

—)■ ^2 o.iJ.{a) G (0, -|-oo). 

a>l a>l 

Then the process converges in distribution in ©(M"*",£^(S')) to the unique solution 

of Smoluchowski’s equation (2.8) started from monodisperse initial conditions co{a,m) = 

/^(®)^{m=l} • 

This justifies part of the title of this paper, in that our model is indeed a microscopic 
model for (the modified) Smoluchowski equation with limited aggregations. 

2.4.2 Precise results in the critical window 

Our second interest is the configuration model. We already saw in the Section 2.3 that 
results on the CM can provide information on our model. More precisely, the divergence 
between our model and a DCM occurs when a large component appears. It therefore 
makes sense to study the largest components in a CM. One naturally expects that, in the 
supercritical case, more than one gelation event will occur, and actually that an order 
a{N)/N will. Hence, we will actually repeat the argument of Section 2.3 over all these 
gelation events, so that we need to be more precise and control the probability with which 
they take place. Moreover, we will have to accurately control the degree of the vertices in 
these components. One of the weakness of the result of [16], Theorem 2.6, is that it is not 
at all quantitative and does not allow this. The strength of this result, however, is that 
the proof can be adapted to show a more precise statement. This is similar to what we 
do for the Erdds-Renyi graph in [24]. 

We shall here consider a CM with parameters n and The important difference 

with Theorem 2.6 is that we do not really want to see these quantities as being n —>■ -|-oo, 
7 r(^) —Tj-j but merely as being part of some “compact” set, and obtain results for fixed n 
and with probabilities which are uniform for all the parameters in this set. 

Let us introduce 

4>n{x) = VkiQiin, 7r("^)))rE^, Wn = Y k^VkiCiin, vr^’^))), 
k>0 k>0 
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and the natural || • ||fc norm on functions on [0,1] or [0,1]^ (which should be clear from 
the context) defined by 


k = Y^ sup D'-f{x) 
2 = 0 


( 2 . 10 ) 


Denote Gn = G^(n) the generating function of its factorial moments, 

and let 

Fn{x) = xG'^{x) = ^ kTr^"'\k)x^, 


k>l 


which is, up to m^\ the generating function of 7r^'^\ Finally, let 


Tj-(n) (n) (n) 

'Jn = J = ^2 ~ 1^1 ■ 


To explain what our “compact” set is, we fix two positive sequences (e„) and (e+) 
with 

<C <1, (2.11) 

as well as three constants r],m,M G (0, +oo). Our main assumptions will be that 

£n ^ In < ^n, "I < ^ TT^"’^(/c)< M. (2-12) 

k>l 


We may then prove the following result. 


Theorem 2.9. Consider a sequence (tt^”)) such that (2.12) holds. Then, for any 5 > 0, 
there exists a constant K > 0, depending only on 6, rj, m, M and (e^) such that, with 
probability greater than 


1 - K 




+ % 


i+v 


we have 


(n) \ / \ 

- 5 \njn< v{ei{n, 7 r^'^'^)) < (2^^ + <5 j njn, ^(62(11, vr^’"))) < 6 njn, 

mn ) \ vrin ) 


and 


f’n 

n7n 



Fn 


<s, 

2 


1 k^Vk{Qi{n, TT^"))) < S. 
^ k>0 


The interpretation of this result is the same as for Theorem 2.6. Let us however list 
the differences with the latter. 

• We require a bit more on jn, namely that it is bounded below by not just 

see (2.11). There is a trade-off here between this assumption and the moments 
assumptions on (vr^"^)). For instance, if we only have £~ 3> the result holds in 

the II • 111 sense, and we have thus an actual strengthening of the result of [16]. On 
the other hand, to get the result in the || • II 2 sense while just assuming e“ 3> 
the proof would require 6-1-7? moments. In any case, since we will anyway need to 
assume (2.2), this result will be sufficient for us. 

• The properties of the components are given for fixed n, not depending on the limiting 
quantities. Specifically, does not need to have a limit, as long as it does not 
become too degenerate. 
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The number of vertices of degree k is controlled very precisely, since this bound in 
the II • II 2 sense translates as 


'^k^ —r;fc(ei(n,7r(’^))) 


-^k7r^'^\k) 


< 5 . 


• The last part of the result allows to control the third moment, but with much less 
precision. It is however enough to say that the measure of degrees of the particles 
which are not in the largest component still have a positive third factorial moment. 

• The probability is uniform for a large set of parameters. 

These considerations will turn out to be of great importance for several reasons. First, 
this result allows us to control precisely which vertices fall and which remain in solution 
after a gelation event. Secondly, we already emphasized that we want to employ this 
result repeatedly: we will use it not only for N, but for all eN < n < N particles (what 
also explain the switch of notation from N to n), as well as for different for which we 
cannot a priori tell whether they converge or not. Finally, this result will be of interest 
to us when 7 „ is of order a{N)/N (see the analysis in Section 2.3 above), so that we will 
use it about N/a{N) times. But in that case, we will get results with probability at least 

I V ^ 1 I 

a{N) \yjNla(N)/NY \ N ) )' 


and this goes to 1 exactly because we assume (2.2). This is the same reason why a stronger 
condition for the threshold {a{N)) appears in [23]. 


2.4.3 Controlling several gelation events: alternative model 

Our final final interest is of course our discrete model. At any given time, the particles 
along with their bounds naturally define a graph structure. Our main goal is to study the 
graph involving the particles in solution. Note indeed that the way that we deal with the 
gel is somewhat arbitrary: it consists of a collection of clusters of size between a{N) and 
2 q;(A^), and we decide that these clusters are inert, but we could also decide to let them 
coalesce among themselves, what would not affect the dynamic. Their behavior is thus 
essentially irrelevant. More importantly, the interesting SOC occurs in solution. 

With large probability. Theorem 2.9 allows to control precisely what happens at the 
first gelation event. The informal idea of the proof is to be able to repeat this argument 
over all gelation events. To make the argument rigorous, we will couple our dynamics with 
that of an alternative model where, at all times, the configuration is that of an appropriate 
configuration model. 

Details are left to Section 7, but in short, the alternative model allows us to see our 
process as pieces of DCM pasted together, every time with different parameters. More 
precisely, we will condition at the configuration in solution after a gelation event, that 
is, Gsink). The first step is to know precisely what this configuration is. To this end, 
denote S{t) the set of particles in solution at time t, and B(t) the total number of 
activated arms in solution at t. For S C [A^j ={1,2,... ,N} and B £ N, we define a 
configuration model CM{S, B) as a uniform choice of B arms of the particles in S, followed 
by a uniform pairing of these arms. Alternatively, we can choose a uniform sequence of B 
arms and link the 2i — 1-th to the 2i-th for i = 1,..., [B/2\. We denote CM'{S, B) the 
graph CM{S, B) conditioned on having no large component, that is we choose a sequence 
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of B arms, uniformly out of all the sequences such that the corresponding pairing does 
not give a large component. Be wary that this is different from choosing first B arms, 
then a uniform ordering out of all those that give no giant component. 

Finally, we say that a stopping time r is a gelation stopping time if, informally, it 
depends only on the particles in the gel. See Section 4.3 for a precise definition. What is 
important is that all the gelation times Tk are gelation stopping times. Our main structure 
result is the following. 

Lemma 2.10. If t is a gelation stopping time, then conditionally on S{t) and B{t), the 
configuration in solution Gs{t) has the distribution of a CM'{S{t),B{t)). 

We can now come back to our reasoning. We want to condition on the configuration 
Gs{Tk) right after a gelation time r^, and recognize this configuration as a DCM on S(Tfc). 
We therefore look at this DCM at the time when we see B{Tk) activated arms, at which 
we exactly have an unconditioned CM{S(vk), B(rk))■ We have then three more steps to 
take. 

1. Check that the conditioning has high probability. This can be done because, thanks 
to Theorem 2.9, we know very precisely the particles that remain in solution after a 
gelation. 

2. Study what happens before the next gelation time To figure this out, the 

previous result allows us to consider our model as a DCM started from an initial 
configuration with the particles S'(rfc). We can then use directly Theorem 2.9, or 
actually the dynamic results of Section 6. 

3. Iterating this reasoning allows to describe our model in between gelation events, and 
will yield the results that we explain in the next section. 

On a side note, there is a similar result to Lemma 4.5 and 4.7, that will not be used in 
the bulk of the proofs, and which gives another description of the structure of our graph. 
Recall that Nt is the number of particles in solution, and tt^ their empirical distribution 
of activated arms. 

Lemma 2.11. Conditionally on Nt and irt, the configuration in solution Gs{t) has the 
distribution of a GM{Nt,TTt), conditioned on having no large cluster. 

In this case, we deal with the classical CM that we explained before. The meaning of 
the conditioning is thus clear: we take a uniform pairing of all the arms, out of all the 
pairings that do not create a large component. We will come back to this result in Section 
2.4.5 below. 


2.4.4 Tightness and limit 


The coupling with the alternative model has two major consequences. First, it allows us 
to precisely control the time between two gelation events. In particular, we will see that 
there exists c > 0 such that, w.h.p.. 


inf T] 
k>0 


(TV) JN) 
k+1 


- r. 


k — 


> C 


a{N) 

1 ^'' 


at least for the in a compact interval. This easily ensures tightness of our processes, 
as we now state. For d > 1, the space is the space of sequences in with finite 

II • ll^-norm, where 

\Mh= + +--- + k^)\u{k,r)\, (2.13) 
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endowed with this norm. With no index || • || = || • ||o. Recall that Pt^\k,r) is the 
concentration of particles that are in solution at time t, with degree r and k activated 
arms. 

Lemma 2.12. Assume that Assumptions 2.2 and 2.3 hold. Then the sequenee 
is tight in and any limiting point is locally Lipschitz-continuous. Moreover, 

for any jd <h + p, the sequence {p^^'^) is tight in D(M+, and any limit point is 

continuous. 


This is the penultimate step towards the description of our system, and it only remains 
to figure out how these quantities {nt) and (pt) evolve. The second consequence of the 
alternative model is that it allows to control the characteristics of the particles in the 
components that fall into the gel — this appeals repeatedly to Theorem 2.9, but with 
different parameters each time. We will find that their generating function solves an odd 
PDE, that is nonetheless fairly easy to solve, see Section 8. The result we end up getting 
is the following. 


Theorem 2.13. Assume that Assumptions 2.2 and 2.3 hold, and denote v 
equation 


(Q(t)-e-*)GUQ(i)) = G,(Q(i)) 


fi. Then the 


(2.14) 


has a unique solution on [Tgep+oo). We denote Q{t) = 1 for t < Tgei, and Q is then 
continuous on 

Moreover, for any (3 < 5 + p, converges in D(M+,M) x D(M+, to a 

limit {n,p) where 

nt = G^{Q{t)) 


and the generating function V’t of pt is 




G,. 


mt) 


- e *)x + e *)y 


We will explain some consequences of this result in the conclusion. Section 9, and 
compute some explicit examples. We will for now be content with explaining how this 
shows SOC in our model. 


2.4.5 Configuration in solution, self-organized criticality 

We now have all the keys in hand to be able to conclude to the graph structure in solution. 
On the one hand. Lemma 2.11 describes this structure. On the other hand. Theorem 2.13 
gives us the parameters of the model. We essentially just have to deal with the conditioning 
in Lemma 2.11, but the alternative model will show that it essentially does not matter. 
The technical details are left to Section 9.1. 

In any case. Theorem 2.13 ensures that —)• tt*, where nt has generating function 


G^rt (x) 


G^ ((Q(t) - e-^)x + e-*) 
~ Gum)) ~ 


(2.15) 


Thanks to Lemma 2.11, we see that to get a sample of our model at time t, we can 
essentially just pick ntiV particles, give them an i.i.d. number of arms according to nt, 
and then create a uniform pairing of the arms. This is a subcritical configuration model 
for t < Rgel) and a critical one for t > Tgel. Check indeed that, by definition of nt, 

miff = (1 — e“*)^m2 < (1 — e~^)mf = 
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for t < Tge\, and that, for t > Tgei, 

= m) - e-^)^G;m)) = (Qit) - e-^)G'^m)) = 

by definition of Q{t). This already proves that our model exhibits some form of SOC. 
However, it is somewhat awkward to deal with CM with random parameters. For instance, 
we may want to say that our configuration is close to a critical CM (with deterministic 
parameters), but this is not really a well-documented notion. 

A nice way around this issue is to use the concept of local convergence that we men¬ 
tioned in Section 2.2.2. We can then give the asymptotic distribution of a typical cluster in 
solution. This is a direct consequence of the alternative model and of the local convergence 
result. Proposition 2.5. 

Theorem 2.14. Assume that Assumptions 2.2 and 2.3 hold. Consider a vertex uniformly 
at random at some time t, and the component Q[yJ{t) rooted at this vertex. Then, in the 
sense of local weak convergence 


Q[^^{t) ^ GW. 




as N ^ -l-oo. 

Since nt is critical or subcritical, these trees are finite, and this could be rephrased by 
saying that, for any finite rooted tree T, 

^hm^P (e(fp)(t) = t) = P(GW,„^^, = T) . (2.16) 

This Galton-Watson tree is subcritical for t < Tgel, and critical for t > Tgel- This shows 
even more clearly that our model exhibits SOC on a microscopic scale. This is similar to 
results in [23]. As a last step, and somehow en passant, we will be able to easily explain 
the intriguing Formula (1.8), answering at last the question of [26]; see Section 9.2. 


3 Smoluchowski’s equation 

The results concerning Equation (2.8) are essentially already proven in the literature, 
though with small differences. This part of the work is different from the remaining of 
the paper, and is merely here to justify that we indeed obtain a model for Smoluchowski’s 
equation. Therefore, we deal with it here once and for all. 

3.1 Well-posedness 

We recall that (the modified) Smoluchowski’s equation (with limited aggregations) is 
written in (2.8), and that Tgei is defined in (2.9). The well-posedness result is Theorem 
2.7. We just need to make clear what we mean by a solution to (2.8). 

Definition 3.1. We say that a family of nonnegative continuous functions {ct{p),p G S) 
is a solution to (2.8) with initial conditions cq G [0,-|-oo)'^ if, for every t > 0, 

(i) 0 < inf^g[o,t](c 5 ,a) < sup^gfo^^j (c*, a) < -hoo; 

(ii) /Q(cs,a)2 ds <-hoo; 
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(iii) for every p G S, 


ct{p) - Ci(0) 



1 


{cs,a) 


(Y^P' ip\p)cs{p)cs{p\p) -YP' PCs{p)Cs{p) 
\p'<p p'SS 


ds. 


Note that conditions (i) and (ii) ensure that the integral in (iii) is well-defined. Up to 
the unsurprising condition (i), this is the definition of a solution given in [26]. Notice also 
that the former has an extra condition in the definition of a solution, but it is unnecessary, 
as was explained for the usual Smoluchowski equation in [23]. 

We will not prove Theorem 2.7. Indeed, due to the similarity between equations (1.4) 
and (2.8), the proof of the existence and uniqueness of the solution is exactly the same as 
in [26], as are Points 1 and 2 of the result. Points 3 and 4 show that both equations are 
just a time-change of one another, and can be checked in a straightforward way. Finally, 
Point 5 is a direct consequence of this time-change. The analysis of [26] also presents a 
representation formula for the solution to (1.4). The same can be done for 2.8. However, 
this gives rise to quite cumbersome formulae, which are not of interest to us. 

It is worth comparing the time Tg^j defined in (1.5), to the time Tgei defined in (2.5). 
We already mentioned that Tggj < -|-oo if and only if Tgei < -|-oo, which is expected as 
Tggj is the gelation time for (1.4), whereas Tgei is the gelation time for the time-changed 
(2.8). Moreover, one can explicitly compute that they are related by the time-change 
described in Theorem 2.7. In our case of interest, i.e. monodisperse initial conditions 
co(a, m) = /i(a)l|^^i|, we also have that mi and m 2 are just the two first factorial 
moments and m^ of p, what makes sense in terms of configuration model, as we 
explained in Section 2.3. 


3.2 Convergence to Smoluchowski’s equation 

Recall that c[^\a, m) is the concentration of clusters with a free arms and mass m in our 
model, that is, there are NcY\a,m) such clusters at time t. Similarly to the result of 
Fournier and Laurengot [14] mentioned above, the only assumption needed for convergence 
is (1.3). The convergence result is Theorem 2.8. 

This result justifies that our model is indeed a discrete model for Smoluchwoski’s 
equation with limited aggregations. Notice also that our model is defined in such a way 
that the cq are monodisperse, but it would be easy to modify it to get any initial conditions, 
and a similar result would hold. However, different initial conditions do not allow us to 
carry out the analysis that will lead to the other results, and which are the main interest 
of this paper. 

We shall now quickly prove Theorem 2.8. The ideas are from [14], already used in [23]. 
There are however two main differences here. 

• There are two parameters a and m, and the coagulations depend on a, whereas the 
threshold is a function of m. This only creates minor issues. 

• Since one has to wait for two arms to be activated to see a coagulation, the process 
(c^'^^) is no more Markovian. The time between two jumps is indeed the sum of two 
independent exponential times and therefore not exponential. The process is then 
called semi-Markov. Once again, this nonetheless allows to write martingales, and 
only minor technical issues occur. For the necessary background on semi-Markov 
processes, see [21]. 
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Proof of Theorem 2.8. Let 


iW _ 


a{N) — l - 1-00 


= E E 

m=l a=l 


aa 


(N) 


(a, m) 


so NA[^'^ is the total number of arms in solution. We always assume in the following 
computations that there are at least two arms in solution (in particular > 2/N), else 
nothing happens. These concentrations evolve as follows. First, one has to wait an expo¬ 
nential time with parameter A^^'^ to see an arm activated. Then, another arm is activated 
after an independent exponential time with parameter A^^'^ — 1, and a coagulation occurs. 
A coagulation occurs between two different p- and pLdusters in solution, with probability 

2ji ■ 

Af\Ar - l/N) 

A p ^ p', and with probability 

P‘PCt [p] 1 P‘Ph \P) 

Ar(AT-l/N)~NAr>iAr-l/N) 

for p = p'. With these probabilities, the new concentrations are , where 

{ \y{p) = \p{p') = -1 if P / y, 

Apy(p) =-2 ifp = y, 

\,P’{P°P') = + 1 - 

The only thing that we are missing is when the two arms are on the same p-cluster, what 
happens with probability 

1 p ■ {p - l)c^P {p) 

where p — k := p — {k,0) for k £ N. In this case, the new concentrations are -|- ;^Ap, 
with 

^p(p) =-1, Ap{p- 2 ) = +l. 

Let E = and /jv = N x {1,..., a{N) — 1}. At any time, the system in then in a 

state rj £ E. For r] £ E, we define 

a{N)-l 

= E E“^(“’"^)' 

m=l a>l 

For any bounded measurable mapping / : E —)■ M, and taking into account symmetries. 
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let us define the generator of (4^^) as 


= 


NA. 


(N) 


+ 


NAij^^ - 1, 


p ■ p'rj(j))p{p') 
P,p'ei%^v (A -1/^j 


E 


f{v + ) - fiv) 


-- E 

N ^ 


p ■ p v{p) 

dN)f AN) 


pTi. - l/AT) 


f{v + - /("?)) 




p- {p- l)v{p) 
— a(^)(aW 1/An 

P&In Ap {Ar, -i/iVj 


/ ( ^ + J^^P ) - /(^) 


4^) - l/(2A^) 


^ p ■ p'p{p)p{p')N (f (p + ^\,pA - fiv) 
\ \ / 


(p,p')e/ 


- ^ E p ■ ?/(p) (^/ (^?? + ^^p.p) - /(^)) 

pSl-n 


+4 • (p -1) (/ (^ + ^^p) - /(^) 


As we mentioned, (c4^) is not a Markov process since the time between two jumps is 
not exponential, but it is however still true [21] that, for every bounded / : i? —)• M, the 
process 

mA'I = nA) - f(4A - f L'«)/{cr>) d. 

Jo 

is a martingale, with quadratic variation 


(M, 


WJ\ _ 


(^(/V) _ 2 /(cf ))l(^)/( cf))) ds. 


Similarly to [23, 14], we can see, by taking 

fiv) = E (« A6)7/(o,m) 


a,m> 1 


for some b > 0, that 


E 


4-) _ i/(2Ar) 


a{N)-l 

E E 


m] 


ds \ <K 


1 4E 

b N 


(3.1) 


for some constant K depending only on the initial conditions. This is the fundamental 
inequality used in the proof, and the rest is just a perusal of the arguments of [23, 14]: 
hrst, tightness is easy to obtain, since there is an order 5N of jumps of size of order 1/N on 
intervals of size <5; then, limits are shown to satisfy (2.8) by using (3.1) and the martingale 
for linear /. □ 
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4 Preliminary results 

This section presents some easy preliminary results. The first is that, for any time, there 
remains a positive fraction of the particles in solution w.h.p., what allows to use asymptotic 
results. The second series of results allow to strengthen convergence in (} to convergence 
in when we have an assumption such as Assumption 2.1 or 2.2. 

4.1 Positive concentration 

Denote Mi the set of probability measures on N, and m[’^^ the subset of /x G Mi such 
that Nfi{k) G N for all A: G N. We consider ^ the set of all couples 

(n, ;u) G [J {n} X M^"^ 

n>l 


such that 

cN < n < N, m 3 > m, ^ < M, (4.1) 

r>0 

and the set of (n, /x) G m ®nch that additionally 


rrio 


m 


^ > 1 + m. 


(4.2) 


We can a priori ensure that the concentration and empirical measure of degrees in solution 
remains in at least on bounded intervals. However, we will also need to show that 

it remains in , and this will require to precisely study how the degrees in solution 

evolve. Typically, assume that half of the particles have degree 3 at time 0, and half have 
degree 2. We could imagine that, we are first going to remove the particles of degree 3. 
When this is done, the empirical measure of degrees in solution will definitely not verify 
(4.2). It turns out that this is not how the system evolves, but this will require more work. 
However, as we mentioned, we have the following result, where we recall that Nt is the 
number of particles in solution at time t, and /xj their empirical distribution of degrees. 


Lemma 4.1. Assume that Assumption 2.2 holds. Then, for any T > 0, there exists 
c,m, M G (0, + 00 ) such that, 

p(vtG[0,r] 


Proof. Since /x^^^ —?■ p, with m 3 > 0, then, by moment assumption, for some 5 > 0 and 
some k >3, p^^\k) > 6 for all N large enough. Denote the indicator function that 

no link on particle i has been activated by time T, and the set of particles with k 

arms. If i G then 

P(Af^^ = !) = (!- e-')^ =: p. 

Since the X^^'^ are independent, by Chebyshev’s inequality. 


P 



Y, Xi-p 
i&lW 



< 


4 

(#/W)V 


E Var(Af)) 

ieiW 


< 


4 

#lWp 
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Since = N > 5N^ then this last term goes to 0 as ^ +oo. Hence, w.h.p., 

at least /2 > p6N/2 := cN particles with k arms have no activated arms, and thus 

are still in solution. Hence, for all t £ [0,r], > c and pt{k) > cN/Nt > c. Therefore 

r>0 * r>0 ^ 

and 

rug* > k{k — l){k — 2)pt{k) > 6c, 

which shows the result. □ 


4.2 Stronger convergence 

The next result ensures that all the weak convergence of measures that we will consider 
can be extended to a far stronger convergence. Recall the definition of and || • ||/3 from 
(2.13). We do the proofs in the case d = 1 for simplicity, but the extension to d > 2 is 
straightforward. 

Lemma 4.2. Assume that is a sequence in such that —)• u in and 

sup < +00 

Ar>l 

fr some 7 > 1. Then, for any ft £ [1,7); ^ u in 

Proof. Let p £ (l,+oo) such that (dp < 7, and let q £ (l,+oo) such that p~^ + q~^ = 1. 
Then Holder’s inequality implies that 


E-' 

r>0 


u^^\r) 


u{r) 


< 




u{r) 



u{r) 


i/q 


and this tends to 0 by the assumptions. 


□ 


This can be extended to convergence of valued functions as follows. 

Lemma 4.3. Assume that is a sequence in such that —>■ u for some 

u £ and that for some 7 > 1 and for all T > 0, 

II Wii ^ I 

sup sup \\ul II 7 <+00. 

Af>lte[0,T] 


Then, for any ft £ [ 1 , 7 ), —)• u in D(M+,£^). 


Proof. It suffices to recall the definition of the Skorokhod topology [7, 19]: for some 
sequence (Xn) of time-changes. 


sup |A7 v(s) — sj -|- sup 

s<T s<T 


(AT) 

u\ , , — Uc 

\n{s) 


0 


for all T > 0. It is then easy to conclude similarly as in Lemma 4.2. 


□ 


Finally, this readily entails that, for processes verifying this 7-th moment assumption, 
we only need to check usual convergence (or tightness) to get a stronger convergence (or 
tightness). 
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Lemma 4.4. Assume that is a process in and that, for all T > 0, there 

is a Mt £ ( 0 , +oo) such that 

P [ sup < Mt] 1 . 

Vte[o,T] / 

Then the following hold. 

1. If —)■ u in distribution in then —?■ u in distribution in 

for any j3 G [ 0 , 7 ). 

2. If is tight in then it is tight in for any (3 G [ 0 , 7 ). 

Proof. The first part follows directly from Lemma 4.3, for instance by using Skorokhod’s 
embedding to assume that —)■ u a.s. The second part means that every subsequence 

of has a subsequence converging in distribution. But then, such a converging sub¬ 

subsequence also converges in distribution in D(M+,£^) by the first part of the result, 
which is exactly what we had to check. □ 


4.3 Combinatorial results 


Recall that S{t) is the set of particles in solution at time t, and that B(t) is the total 
number of activated arms at this time. A configuration model CM{S,B) is defined by 
picking a uniform sequence of B arms of the particles in S and joining the first and the 
second, the third and the fourth, and so on. The conditioned graph CM'{S,B) is this 
graph conditioned on having no large component. Precisely, it is given by the pairing 
corresponding to a sequence of B arms in S, chosen uniformly among all the sequences 
such that the corresponding pairing creates no large component. The main combinatorial 
result is the following. 

Lemma 4.5. For any t>Q, conditionally on S{t) and B{t), the configuration in solution 
Gs{t) has the distribution of a CM'{S{t),B{t)). 

Proof. Denote by a{t) the (ordered) sequence of arms activated up to time t (whether 
they are in solution or not), and |cr(t)| its length, that is, the total number of activated 
arms. Define S{a) to be the particles in solution when we activate the arms in the order 
a and perform our algorithm, that is, bind the first and second arm in solution, the third 
and fourth in solution, and so on. For S C [N], let also crs{t) be the sequence of arms 
activated on particles of S. 

For S C [N] and A: > 0, denote Ls^k the ordered sequences £ = (£ 1 ,..., £fc) of fc distinct 
arms in S, and Lgf. those sequences such that linking £i and £2, £3 and £4 and so on, does 
not create a large component. Our goal is exactly to show that 


F{S{t) = S,as{t) = £) 

is independent of ^ G ^. 

To begin with, we can split this event according to the arms that are activated in 
S = [A']\5', by writing 


P(5(t) = .S,CT5(t)=A) 

= E E ^iS{t) = S, as{t) = £, a^it) = £) 

= E E ^{W{t)\ = k + r)F(^S{t) = S,as{t) 


£, a-o{t) = £ |cr(t)| 
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Now, as we already mentioned, conditionally on \cr{t)\, the configuration is given by a 
uniform ordering of |(T(t)| arms, i.e. a{t) is uniform in L[7v],|CT(t)|- Therefore, for I G Lg 


P ( 5 (t) = S, as{t) = i, ag{t) 


a{t)\ = k + r^ 

_ #{ 0 - £ T[jv]^fc+^,crg = i,ag = I,S{a) = S} 

i^{(T G I-'[N],k+ri — ^1 


The denominator clearly only depends on k and r. Therefore, it suffices to show that the 
numerator is independent of So consider another G L'sk- There is a clear bijection 


cj G L 


[N],k+r 


,as = £,ag = £,S{a) = |cr G L[jv],fc+r, 0-5 = = ^,S{cr) = . 


Indeed, take an element a of the former set, and replace ii by for i = 1,... ,k. This 
obviously provides a u' such that a'g = i',ak = £, and minute of reflection also allows to 
check that S{a') = S. There is an obvious inverse mapping, and we have thus a bijection 
between these two sets, what allows to conclude. □ 


As in [23], it is possible to prove an extension of this result. We first define a gelation 
stopping time as in [23]. For any t >0, the natural filtration of our model is 

Tt = CT a S ^ ) 

where A is the set of arms and Ca the clock on arm a. Similarly, for a subset S of [A^], we 
define the filtration generated by clocks attached to particles in S 

3'f := cr (eal{e^<t},a G A 5 ) , 

where As are the clocks on the arms of the particles of S. 

Definition 4.6. We say that r is a gelation stopping time if 

• r is a (Ti)t>o-stopping time, 

• for any S C [N] and f > 0, conditionally on {S{t) = S}, rl{.r<t} is independent of 

3-f- 

Two important examples of a gelation stopping time are 

• any given deterministic time t > 0. 

• the fc-th gelation time for a given A: G N. 

The fact that the gelation times Tk are gelation stopping times is easy to see. In¬ 
deed, conditionally on {S{t) = S}, is determined by the k-th. gelation time in 

(Ss(s))s<4) and it is therefore independent of . 

On the other hand, for instance, the first time after that a component of size at 
least a{N)/2 is created is a (Ttj-stopping time, but is not a gelation stopping time. 

Similarly to what is done in [23], the previous result extends to gelation stopping times. 
The proof is similar to that of [23]. 

Lemma 4.7. If t is a gelation stopping time, then conditionally on S{t) and B{t), the 
configuration in solution Gs{t) has the distribution of a CM'{S{t),B{t)). 
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Instead of considering the particles in solution themselves, we can take a lighter con¬ 
ditioning on the number of particles in solution and their degree distribution. We thus 
define a CM{M, B) by taking M vertices with degree distribution // and then taking 
the pairing given by a uniform choice of a sequence of B arms. As above, CM’{M, fj,, B) 
is a CM{M, fi, B) conditioned on having no large component. By summing the previous 
result over all the S with the same cardinality and the same measure of degrees, and by 
exchangeability, we get the following. 

Lemma 4.8. If t is a gelation stopping time, then eonditionally on and pL^, the 
eonfiguration in solution Gs{t) has the distribution of a CM\Nr, Ht, 

Conditioning on the total number of activated arms in solution will turn out to be 
the most useful to us; see in particular the description of the alternative model in Section 
7. However, we can also condition on vr^, the empirical distribution of activated arms in 
solution, to get a more natural result directly related to the CM. The proof is again done 
in a similar fashion. 

Lemma 4.9. If t is a gelation stopping time, then eonditionally on Nj- and tTt-, the 
eonfiguration in solution Gs{t) has the distribution of a GM^NrjTTr), conditioned on 
having no large cluster. 

This more natural result will only be useful to us when we want to study typical 
clusters in solution and use Proposition 2.5, see Section 9.1. 

5 Largest component of a slightly supercritical 

CM 

This long section is devoted to proving Theorem 2.9, and providing some extensions to 
vertices with “types”, which should be thought of as the number of free arms of our 
particles. This section is purely about random graphs, and does not use any specifics of 
our model. We recall that the notation and the statement of the main result are in Section 

2 . 2 . 

5.1 Proof of Theorem 2.9 

We outline the argument of [16] for proving Theorem 2.9, and point where it can be slightly 
improved. The beautiful idea in [16] is to perform simultaneously the random uniform 
pairings of half-edges and the exploration process of the components, in the following 
way. Color^ the vertices as white or red and the half-edges as white, yellow or black; 
white and yellow half-edges are also called bright. Start with all vertices and half-edges 
white. Pairing one edge with another will be decided by giving the half-edges i.i.d. random 
maximal “lifetimes” Tx with distribution £(1): namely, each half-edge spontaneously turns 
black with rate 1 (unless it was colored black earlier). 

Cl If there are yellow half-edges, go to C2. If there are no yellow half-edges (as in 
the beginning), select a white vertex by choosing it uniformly at random^, color it 

^In their article, Janson and Luczak rather label vertices as sleeping and awake, and half-edges as sleeping, 
active and dead. We use colors here as to not confuse the exploration with the dynamics of our model. 

'^Janson and Luczak rather choose a white half-edge uniformly, but this only changes the order in which we 
build the components: size-biased here, biased by their number of half-edges in [16]. This slight modification 
has the advantage of simplifying a bit the end of the argument. 
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red and all its half-edges yellow. If the chosen vertex has no half-edges, repeat Cl, 
otherwise go to C2. 

C2 Pick a yellow half-edge e (which one does not matter) and color it black. 

C3 Wait until the next half-edge e' turns black (spontaneously), and join it with e to 
form an edge of the graph. If the vertex v to which e' belongs is white, change v to 
red and all other half-edges from v (but e') yellow. Repeat from Cl. 

Let W{t), Y{t), B{t) = W{t) -|- Y{t) be the number of white, yellow, and bright 
half-edges at time t. We naturally define these processes as cadlag. 

As observed in [16], this algorithm allows to construct the configuration model compo¬ 
nent-wise, and components are created each time Cl is performed. Note also that the times 
Cl is performed exactly corresponds to the times where Y cancels. Thus, the exploration 
of the successive components are performed exactly during the successive excursions of Y 
above the origin. 

Remark 5.1. On the other hand, given a graph 9, this algorithm also allows us to explore 
9 eomponent by eomponent. Consider indeed that eaeh edge is made up of two half-edges, 
and perform the same algorithm, up to two slight differenees: forget about the exponential 
cloeks, and replaee C3 by 

C3’ Take e' the other half-edge attaehed to e. If the vertex to whieh ^ belongs is white, 
ehange this vertex to red and all other half-edges from this vertex (but e') yellow. 
Repeat from Cl. 

It is elear that, if 9 is random and obtained by the eonfiguration model, then the sequenee 
of components explored by this algorithm has the same distribution as the sequenee of 
components constructed by the original algorithm. Since each exploration of a new com¬ 
ponent is started by choosing a vertex uniformly at random, then the order in which we 
build the components is actually biased by their size. 

Let us fix to > 0. In the whole proof, C denotes a constant depending only on to and 
on the fixed parameters S, M, m,e^, and which may change from line to line. 


1. Observe that {B{t)) is a death-process starting from — 1. B{t) is decreased 

by 2 at rate B{t) (except if there remains only one bright half-edge). Lemma 6.2 in [16] 
ensures that 


sup 

B{f) — exp(—2f) 

2 

7<7ni0 




< Cnjn 


(5.1) 


Now, let Vfc(t) be the number of vertices with degree k > 1 having all their half-edges 
with lifetimes Tx > t. Again it is a death process, decreasing by 1 at rate kVj.(t), starting 
from mr^"'\k) (or mT^'^\k) — 1 for one of them), and as in (6.5) of [16], we have 


E 


sup 

t<7ni0 


Vk{t) — mr^'^\k) exp(—fct) 


< Ckn'~fn'^^'^\k), A: < l/7„ 


(5.2) 


The interest in 'Vfc(t) is that, at least at the beginning of the construction, it is a very 
decent approximation of 14 (t), the number of vertices with degree k > 1 that are still 
white at time t. Better yet, W(t) = J2k>o^^kit) is also a very decent approximation for 
the number of white half-edges W{t) = J2k>o again at least at the beginning of the 

exploration. More precisely, one should observe that the difference can only come from 
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half-edges which have been colored at the first step Cl of the algorithm, therefore W and 
W should remain close as long as we have not explored too many components. 

Slightly precising the argument of the proof of Lemma 6.3 in [16], we find that 


E 


sup 


W(t) — n ^ k'K^'^\k) exp{—kt) 


k>l 


l/7n 





- 

< y A:E 

sup 

Vk{t) — n7r^"'\k) exp{—kt) 

+ T, kE 

sup 

Vk{t) — n7r^'^\k) exp(—fet) 

k=l 

t'^'y-nto 


k>l/'yn 

_t<'y„to 

. 


. - 

— ky Ckn'yn'^^'^\k) + ^ km^^'^\k) 

k=l k>l/'yri 


1 /7»i ! - 

< EV k^'K^^){k)-\-n ^ kTr^^\k) 

k=l k>l/'y„ 


< 


/l/7n 

CV^ E k^+^7r(^\k) 

\k=l 


1/2 /l/^„ \ 1/2 

+117^" E k^^^TT^^Kk), 

,k=l ) fe>l/7„ 


where, at the second line, we use (5.2) and Cauchy-Schwarz for k < l/7n, as well as the 
trivial bound 'Vfc(t) < mT^'^\k) for k > l/jn- Finally, our moments assumptions (2.12) 
allow to conclude that 


E 

sup 

'W(t) — n y k'iT^'^\k) exp(—/ct) 


i<7nt0 

fc>i 


In fact, from the proof, it even holds that 


< C'(Vn> + ?^7n’^’i)- 


(5.3) 


E 

y k sup 

'Vfc(t) — mT^'^\k) exp(—/ct) 


fc>l t<7nt0 



which will turn out to be useful later on. 




(5.4) 


2. Let now y(t) = B{t) — W(t), which should be a good approximation to Y{t) = 
B{t) — W{t), at least at the beginning of the exploration, and define /i„(t) = — 

J2k>i k'^^"'\k)e~^^. Then (5.1) and (5.3) readily imply that 


E 


sup 

t<tQ 


-7n^y(7n/) 

n 


< C 




+ 7^^ 


(5.5) 


Note that hn{t) = F„(l)e 2*—ir„(e *). Since 7n < 0) compute the following, 

where the constants hidden in the O(-) depend only on (e^) and M, and are uniform in 
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t G [0, to]: 

hnilnt) = 

= F„(l)(l - 27 „t + 2^lf + 0{jI)) - [f„( 1) + F;(l)(e-^"‘ - 1) 

+ ^F"(l)(e-^"* - 1)2 + 0(73)- 

= -^n(l)(l — “^Int + + 0(7n)) “ -^n(l) + -^n(l) ^“7ni + 

+ “*■ ^(7n))^ + 0 ( 7 ^) 

= 7ni (|— 

+ Umi^ - ^ ^ + 0 ( 7 ]^) 


= 7ni + 7ni^ 

= 7ni (^1 - +0(7n)> 


where we use the definition of the and 7 ^ = — m^\ In short, we get 


In'^Kilnt) -t(^l- 
Note that (5.6) is a purely deterministic statement. 


sup 

t<to 


< C'7„ 


(5.6) 


3. Now, observe that I4(t) < 'Vfc(t), so that W{t) < W(t), and that moreover, W — hF 

can only increase when Cl occurs. Take such a time t when some half-edges are colored 
due to Cl occurring. Then Y{t~) = 0 and Y{t) is the number of vertices of the vertex 
chosen at Cl, so that Y{t) < max{i, 7 r("')(i) 7 ^ 0} < (Mn)^/^^'''^) < Cn.'j^, where the first 
inequality comes from (2.12), and the second from the fact that e~ We deduce 

that 

W{t) - W{t) = W(t) - B{t) + Y{t) < W{t) - B{t) + Cn-fl 
At other times, W — IT decreases, thus 

W(t) — W{t) < sup(W(t) — B{t)) -I- Cwjn = — inf y(s) -|- Crt'y^. 

s<t 

For the same reason, up to another factor max{z, 7 r(”^(z) 7 ^ 0} < Cnj^, this even actually 
holds until T*', the first time strictly after t when Cl is performed. To conclude, for all 
t > 0, 

0 < sup (W(s) — IT(s)) = sup (T(s) — y(s)) < — inf V(s) -|- Cn'y^. (5-7) 

s<t+ s<t+ 

4. Fix now d > 0 and 5' = 5/{2M). Define tn = 2jm^\ and notice that Pn{t) = 

t(l — /2) is positive on (0,t„) and vanishes on the boundary. For e > 0 small 

enough, it is clear that the parabola Pn may be included in a tube of height e as in Figure 
1. In other words, with ± e, P~ vanishes at two points, one on [0, YI2\ and one 

on [tn — 5'/2,tn], whereas P^ > 0 on [0,tri] and F+(t„ -|- 6'/2) < 0. Moreover, one can 
choose £ uniform as long as tn is in a compact interval of ( 0 , -|-oo), what holds from our 
assumptions. 
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Figure 1: The parabolas , Pn, and P^. 


Consider the event 


Er, 



7n hni^nt) 



From (5.5) and Markov’s inequality, we have IP(F„) > 1 — Pm where 


4C / 1 

Pn = — [ 


1+r] 


From now on, we work on E^, and take n large enough (which can be chosen uniformly) 
such that C^n < e/4 and < e/4. We thus have 


in/ -7n y(7ns) > inf -7n Kilns) - - > inf P„(s) - C-fn - - = -Cjn - T > --• 

S^tn Tt S^tn Tl 4 S^tfi 4 4 


where the second inequality stems from (5.6). With (5.7), this implies 


sup -7n^ l^(7ns) - y(7nS)| < “ in/ -7n^y(7nS) + C^in < 
n s<tn n 4 


(5.8) 


and from (5.6) and since we work on En, we get 

sup -7“^y(7nS) - Pn{s) 
s<t+ ” 


3e „ q 
< ^ + Cn-f^ < e. 


(5.9) 


In other words, on [0,t/^], s e-)- n“^ 7 “^y( 7 „s) is stuck between P~ and P^. Hence, there 
must be an excursion of y( 7 „-) above 0 on an interval [ri,T 2 ] containing the interval 
[y/2, tn — 5'/2]. We have Ti G [0, d'/2]^ and T 2 is the first time after tn — d'/2 when Cl is 
performed. If T 2 > tn + d'/2, then + d'/2, so n~^^~‘^Y{'yns) is below P+ at least 

on the interval [0,tri + d' / 2], This forces Y{tn + d' / 2) < 0, which cannot happen. Hence 
T 2 G [tn — d'/2,tn+d'/2]^ and there is an excursion of y with size in [7n(in —7n(^n+'^0]- 
This corresponds to a component Cq. 
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5. Take now jnTi < 7 nT 2 < 'Jntn be two times when Cl is performed, so that the 
exploration of a component starts at 7n'7‘i, then one or several components are explored, 
and an exploration ends at 'yn'T^■ Let us compute the size and number Wk{Ti,T 2 ) of 
vertices of degree k in these components. We will first compute the number of vertices 
Wkiji, T 2 ) in these components. 

To begin with, define 

N{t) = j2Vk{t), m = 

k>l k>l 


Recall that I 4 < so that T — 'y = W — 1T = X) k\'^k ~ Lfc|, and (5.8) then implies 


sup -7„^ |iV(7nt) - W(7ni)| < 
t<tt ^ 4 


Moreover, (5.4) and Markov’s inequality yield 


sup 

t<t+ 


-In - 7n ‘^9n{e '^"*) 

n 


- 4 


(5.10) 


(5.11) 


with probability greater than 1—pn- Clearly, Wk{Ti,T 2 ) is the number of vertices of degree 
k which have turned from white to red during [ 7 nTi, 7 nT 2 ), i.e. Wk{Ti,T 2 ) = 14 ( 7 n'ri) — 
VkilnT^), so that N{ti,T 2 ) = A^(7nn) - N{'^nT 2 )- Since ti,T 2 < (5.10) and (5.11) 

imply 

\N{n,T2) - (5'n(e“'^"^0 “ 5n(e"^"'^2)) | < 2£. 

A first order Taylor expansion then yields 

ffn(e“^"^') - ffn(e“'^""'") = 7n5r(,(l)(T2 - Ti) + 0{'yl) = Jnm^^\T2 - Ti) + 


where O(-) is uniform in the parameters. The last two equations then yield 


n 


7 n N{ti,T 2) - iT2-Ti)m\ 


(n) 


< C'^n 


We may from now on consider n large enough (uniformly in the parameters), such that 
C'fn < (^/2. Therefore, we finally obtain 


-In ^N{ti,T2) - {t 2 - Ti)mi 
n 


(n) 


< S/2. 


(5.12) 


In particular, for ri = Ti and T 2 = T 2 , we get, remembering that T 2 — T 1 G [tn — S',tn + (5'], 
that the number of vertices N{Ti,T 2 ) of |Co| verifies 


1 

n 


7^'|eo|-2 


m 


(n) 


m 


(n) 


< 5/2 + MS' < S. 


(5.13) 


6. We have therefore found a component Cq of the right size. We shall now prove that 
with probability at least 1 — Cpn there is no larger component in the graph. 

First, we know from the above reasoning that, with probability at least 1 — 2p„, the 
following happen. 

• We build, on the interval [Ti,T 2 ), a component Cq of size in {sn — S,Sn + S)n'yn, 
where Sn = 2rn\'^^/rrL'^\ 
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Any component built before Co has size less than 


n'yn < nyn = dn-fn- 


These equations come from (5.12) and (5.13). We shall call this event £i. Consider further 
the events 

• 8 , 2 - we build a component of size {sn — 6,Sn + S)n'jn and at least another one of size 
greater than Sn'jn] 

• £ 2 - ''^6 build at least one component of size greater than Sn'jn before a component 
of size {sn - S,Sn + 5)n7„. 

Obviously £'2 C £2, and £'2 n £1 = 0 . 

As mentioned in Remark 5.1, the components are built in a size-biased order, and thus 

IP(£2|£2) > S/{Sn + 2S)>p 


for some p > 0 uniform in the parameters. Hence, 


ip(£2) 


p(£' 2 n £2) 

P(£'2|£2) 


< ^ip(£ 2) < ^(1 -iP(£i)) < 

p p 


2 pn 

P 


Finally 


T(£i\£ 2) > nSi) - ^(£ 2 ) > 1 - 2p„ - ^ = 1 - 2 f- + 1 ) p„. 

p \P J 

The first part of the result follows after noticing that £i \£2 is exactly the event considered 
there. 


7. Let us write w.h.p. to mean with probability at least 1 — Cpn, for some C uniform in 
the parameters. The reasoning above shows that, w.h.p., Cq is indeed Ci(n, vr^^^). Hence, 
to complete the proof, we need to show that 




1 


n-Tn 


Vk{Qei) 


-^k7r^'^\k) 


<5, 


- k^VkiQo) < b 


(5.14) 


w.h.p. 

To begin with, w.h.p., Cq is the largest component in the graph, and has size at least 
cw-fn, where c > 0 is uniform in the parameters. Since the components are constructed 
in a size-biased order, there is probability at least cny^/n = cjn that it is explored at 
the first step. If not, there is probability at least cn'^n/n' > c^n that it is explored at the 
second step, where n — n' \s the size of the first explored component. Clearly, this shows 
that the number K of times that Cl is performed before exploring Cq is stochastically 
dominated by a geometric distribution S(1 — C 7 n) with success parameter 1 — 07 , 1 . So take 
[5 > 1 such that 2 + [5 > 2/[A + p). We have then 

P (g(l - C7n) > l/7n’^^) = (1 - < exp -c/p^ < . 

Therefore, w.h.p., K < 7n^^"*'^^. 
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Now, recall from Point 3 above that 0 < Vfc(t) — 14 (t), and that this can only increase 
when Cl occurs. Recall that So is explored on an interval [7n7i,7n72)) with T 2 uniformly 
bounded. Then, w.h.p.. 


sup < (max{z,7r(’^)(z) / 0 }fK 

< 2 '^ln, 


for n large enough, where the last inequality comes from the fact that 'Jn ^ ^ 

and the choice of (3. This is when we most crucially use this assumption. 

A similar computation as in Point 1 of the proof allows to obtain 


y 

sup 

Vk{t) — nA"'\k) exp(—fct) 

k>l 

_t<7„T2 


Markov’s inequality therefore ensures that 

p I y /c^E 

sup 

Vk{t) — nA"'\k) exp{—kt) 

\fc>i 

t<1nT2 

. 


<C[J— + nj, 


> 2™7n 1 < C 


,2+r) 

n 


(5.16) 




+ 7 


l+r) 

n 


Along with (5.15), this gives 


sup ^ Vk{t) — mr^"‘\k) exp{—kt) < drijn 


w.h.p. But i;fc(eo) = VkilnTi) - Vk{lnT 2 ), so that 
y k"^ —Ufc(eo) - —TT^^\k) 

- n7n 0'^ V / 


A:>1 


In 


< 26. 


A simple Taylor expansion then allows to conclude to the first part of (5.14). To get the 
second part of (5.14), note that similar computations allow to obtain 


y k^ -VkiQo) - A^\k) 


fc>i 


< 26, 


so that 


y k^Vk{Qo) < 2(5 + y k\^^\k) 


k>l 


k>l 


< 2(5 + C 7 n y A:V”)(A:) 


k>l 


and the assumptions (2.12) allows to conclude. 

All these reasonings are for n large enough, but this can always be done uniformly in 
the parameters. It therefore suffices to take a larger K in the statement to deal with the 
small re. 


Remark 5.2. A estimation of 4>n could be direetly obtained from Point 5 of the proof. 
However, to get bounds on the seeond and third moments, we need to use the slightly 
different teehniques of Point 1, essentially due to the faet that the techniques leading to 
(5.10) care only provide a bound for the first moment. 
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5.2 Extension to particles with free arms 

5.2.1 Result 

In our model, coagulations occur thanks to the activated arms, but each particle also 
carries a number of free (i.e. not activated yet) arms. Hence, what happens at a gelation 
event only depends on the distribution of activated arms. However, given the state in 
solution after a gelation event, the dynamics depend on the number of free arms. For this 
reason, we will need to keep track of the number of free arms. More precisely, we will 
follow the number of particles with a certain number of active and free arms, and it makes 
sense to first study this for the configuration model. Rather than the number of free arms, 
we will instead use the degree as the parameter, since it is a constant for a given particle. 
Vertices with k activated arms and degree r will just be called (/c, r)-vertices. 

Therefore, let us consider n vertices, with (/c, r) of them being (/c, r)-vertices. Only 
the active arms are involved in the coagulations, and we therefore define as 

a uniform pairing of the active arms, and let and be its largest and 

second largest component. For k,r > 0, let 

vrW(fe) = J2&Kk,r), /iW(r) = -j2&\k,r). 

r>0 k=0 


Therefore, if one forgets about the free arms, then has exactly the same 

distribution as a as in the previous section. We further define 


y) = XI ^ r)x^y\ 

r>0 k=0 


^n{x,y) 


= X 


d&n 

dx 


{x,y), 


as well as = m, 


-(") 


t{'») 


and 771 = 7^^ • We reserve this notation with Gothic letters for 

the “effective degree”. This will become more relevant for Theorem 5.6 below. 

Under the same assumptions on {ttu), Theorem 2.9 is then still valid for 
However, it does not tell us about the distribution of the degree of the vertices that 
belong to the largest component. For this reason, let us denote by the number of 

{k, r)-vertices in a component C. Introduce 


'ifnix, y) = X X 
r>0 k=0 


X rVr(ei(n,^W)). 

k,r>0 


We shall need assumptions corresponding to (2.12). Note that the “4 + rj moment” as¬ 
sumption on is implied by a similar assumption on However, we will need to 

assume 5 + t] moments here. We shall therefore fix two positive sequences 

<C <C 1 


and three constants r],m,M G (0,-|-oo), and assume the following. 

Assumption 5.3. 

, fn < X < M. (5.17) 

k>l 

The extension of Theorem 2.9 corresponding to this model is the following. 
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Theorem 5.4. Consider a sequence such that (5.17) holds. Then, for any 5 > 0, 

there exists a constant K > 0, depending only on 6, rj, m, M and (e^) such that, with 
probability greater than 


1 - K 






we have 


and 


.nr 


(n) 


m 


(n) 


- n7n < |ei(n,^(’^))| < ^2^^ + (5^n7„, |e2(n, < (5n7„ 


- fjn 

n-^n 



S'n 


2 


-Wn < S. 
n 


Of course, the only difference with the previous result is the last part, concerning the 
number of (fc, r)-vertices in the largest component. It says as before that the vertices that 
make up this largest component are chosen according to a size-bias of their number of 
active arms, whereas the degree is irrelevant. 

In particular, as the free arms do not matter for the coagulations, this result could be 
rephrased using abstract types rather than the total number of arms, but at the price of 
slightly awkward assumptions. Specifically, we systematically use that k < r. 

A final note is that, because of the 5 - 1-7 moment assumption, rather than 4 -|- 7 in 
Theorem 2.9, the lower bound on 7 ^ could be slightly improved. For all intents and 
purposes, this result will be more than sufficient. 


5.2.2 Proof 


The proof relies on the exact same argument as the one of Theorem 2.9. So we perform 
the algorithm Cl, C 2 , C3 on the active arms, and consider all the same variables and 
notations. Note that the assumptions (2.12) on are in force, and thus every single 
step of the proof is still valid. We wish to prove that 


EE’ 

r>0 fc=0 


n'yn 


-Vk: 


{eo)-^knC^^\k,r) 


m. 


< d 


(5.18) 


w.h.p. 

To this end, we need to keep track of the extra information provided by the total 
number of arms, so we additionally introduce the variable Vk^r{t) for the number of {k, r)- 
vertices that have all active half-edges with lifetimes Tx > t, and similarly Vk^r{t) for the 
number of {k, r)-vertices that are still white at time t. 

The only difference comes from Step 7 of the proof of Theorem 2.9. First, 0 < — 

Vk^r: this can only increase when Cl occurs, which happens K <^n times w.h.p. 

Moreover, w.h.p., the largest component Co is explored on an interval [ 7 nTi, 7 nT 2 ) with 
T 2 uniformly bounded, and thus 

sup r'^\Vk,r{t) - Vfc,r(i)| < (max{ 7 //fo)(z) / 0})‘^K 

< dn^ri- 


This is just (5.15) in this context. 
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We now need an equivalent of (5.16), which we shall get as in Point 1 of the previous 
proof. As for (5.2), Vk^r is a death process decreasing by one at rate kV^^rif) and starting 
from We therefore have 


E 


sup 

t<'rnT 2 


r) ex.Y>{-kt) 


<Ckn'yn&'\k,r), A: < l/y^. (5.20) 


sup 

VkAt)-n^^^\k,r)e-’^^ 

t<~fnT2 



sup 

t<lnT 2 


Vfc,r(i) - nC^'^\k,r)e 


— kt 


sup 

VkAt)-A^^\k,r)e-’^^ 

_t<7„T2 



As before, we write 

E 

A;,r>0 

l/7n r 

< 7n^ E ^E®^ 

r=0 fc=0 

+ E ^^E®^ 

»'>l/7n k = 0 

and deal with each sum on the RHS similarly. First, the second sum on the RHS is 
bounded by 

r 

r‘^'^n^^^\k,r) = n ^ r^//^"^(r) < ^ ■ 

r>l/^n k=0 r>ll')„ r>l/')n 

As for the first sum on the RHS, we use (5.20) and Cauchy-Schwarz to write that it is 
bounded by 

l/7n r _ l/7n / r \ 1/2 / r \ 1/2 

Cy/n^ E E V !■) < C^/r^ E ^ E ^ ) E ) 

r=0 fc=l r=0 \fc=l / Vfe=l / 

l/7n 


< Cyfr^ ^ 

r=0 

/l/7r> \ /l/7»i 

< E E 


1/2 


, r=0 


r=l 


< C^/n^, 

where the last step comes from (5.17). To summarize, we obtain 


E 


7 sup Vk,r{t) - n^^'^\k,r)e 

k,r>0 l< 7 nT 2 




(5.21) 


This is just the equivalent of (5.16) in this context. Formula (5.18) and the rest of the 
statement can then be obtained as in the previous proof. 


5.3 Percolation on the configuration model 

5.3.1 Result 

In our real model, at each time t, arms are activated independently with probability 
Pt = 1 — e“b Hence, the degree of a particle is fixed, but its number of activated arms is 
random. We will thus study the corresponding configuration model in this context. 
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To be precise, assume that there are n vertices, of them being of degree r. 

Let Gn = G^(n) be the generating function of and define = mf the factorial 
moments of Let pn G [0,1] and assume that each arm is activated independently 

with probability pn- We then independently perform a uniform pairing of the activated 
arms, to get a configuration model GM{n, p,^'^\pn)- We naturally denote Qi{n, p^'^\pn) 
and Q 2 {n, p!''^\pn) its largest and second largest component, Vk,r{G) the number of {k,r)- 
vertices in a component C, and finally 

r 

'^n{x,y) =^^Vk,r{Ql{n,p^'^\pn))x^y'", Wn= (^1 ' 

r>0/c=0 k,r>0 

Of course, conditionally on the activated arms (in a “quenched” setting). Theorem 5.4 
will apply. Two issues however remain: first, one needs to check the assumptions (5.17). 
This can be achieved, as we will see, when pn is chosen so that we get a critical model, 
namely 

(n) 

Secondly the result will then provide information that is conditional on the activated arms, 
that would not be convenient to use. However, by the law of large number, the initial 
distribution should be almost deterministic. Indeed, if we let 


[p]1 =[Ap\^-p) 


r—k 


(5.22) 


then should be about y^'^\r)[pnYk- Therefore, the generating function of 

^(n)(^,^) giiould be close to its “annealed” version 

&n{x,y) = f^^"^\^)[Pn]kx’'y'' = GniiPnX + 1 -Pn)y)- 

r>0 fc=0 

Ignoring the first order terms, it makes sense to define 

(rn{^ 

qn{x,y)= 1 -^x + l-^\y. 


, rrio 


rUn 


With a slight abuse of notation, we will write x&n for the function x e-)- x&n{x), and so 
on. We shall again fix two positive sequences 


n <C <C 1 


(5.23) 


and three constants r],m,M G (0,+oo), and do the following assumptions. 

Assumption 5.5. 

m < < mY\ m < < M. (5.24) 


k>l 


We will prove the following consequence of Theorem 5.4. 
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Theorem 5.6. Consider sequences and {pn) such that (5.24) holds. Then, for any 

5 > 0, there exists a constant K > 0, depending only on 6, rj, m, M and (e^) such that, 
with probability greater than 


l-K 



we have 


and 


and 




m] m. 


m 


(n) 


-d]nen< |ei(n,//(”\pn)| < ( 2 




m] m, 


m 


(n) 


+ S nSr. 


\e 2 {n, Pn)\ < 5ne, 


- ifn 

nSri 


(n) 

' (n) 

run 


xgG'Jqn) 


<S, 

2 


-Wn < S. 
n 


This shows that, in a configuration model with parameters n and the right per¬ 
colation threshold to obtain a largest component of size of order nSn is rn'i'^+ en). 
Moreover, taking x = 1 or y = 1 in the second part of the statement shows that the 
particles belonging to the largest component are chosen both by a size-bias of their num¬ 
ber of activated arms, and by a size-bias of their total number of arms. This will prove 
fundamental in our analysis. 


5.3.2 Proof 

Let us proceed with the proof of this result. As we mentioned, the number nf^^'^\k,r) of 
(A;,r)-vertices should be close to p^^\'r)[pn]k- More precisely, we will prove the following 
inequality. We denote here for the empirical distribution of {k, r)-vertices when each 
arm is activated independently with probability p, and Q for the probability under which 
the activated arms are chosen. 


Lemma 5.7. For any 5 > 0, there is a constant K G (0,-|-oo), depending only on 5 and 
M, such that for any p G [0,1], n > 1 and > 0, 


EE"' 4")(fe,r) - p^^\r)[p\l >5an\<K[ + 


\.r>0 k=0 




nay 


(5.25) 


and 


f^;^\k,r) - p^-\r)[p]l >5 <K[aY^ + ^^ 


(5.26) 


We delay the proof of this result to the end of this section, and shall first use it to 
prove Theorem 5.6. 


Proof of Theorem 5.6. Let us use the same notation 7n, ®n, Sn as in Theorem 

5.4. Note that these quantities are random, but conditionally on their values. Theorem 
5.4 applies, provided the assumptions (5.17) hold. We shall therefore start by checking 
them. 
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To begin with, the last assumption of (5.17) is just the same as in (5.24). Now, from 
(5.25) and (5.26) with an = Sn, p = Pn and 6 small enough (to be fixed later), we have 


&Hk,r) - p(^\r)[pn]l 


r>0 k=0 


< 6 er, 


k,r>0 


&\k,r)-p(^\r)\pn]l <6 (5.27) 


with probability at least 



As usual, we will just write w.h.p. for such a probability, where K depends only on 
the fixed parameters, and C will be a constant which may change from line to line, but 
depends only on the fixed parameters (but not 5). 

To begin with, the first equation of (5.27) implies that 


^>1 k=l 


is, up to 5£n, 

E E kp^''\r)[pn]l = '^Pnrp^'^\r) = pnmf\ 

r>lk=l r>l 


Similar computations hold for the second and third factorial moments, and thus 


(n) (n) 

m\ - pnm\ 


< C5enj 


in) 2 M 

- Pnm-2 


< C 6 £, 


(n) 3 (n) 

- Pn^s 


< C5. (5.28) 


A couple easy Taylor expansions then ensure that the corresponding 7 „ verifies 


7n 


m, 


(n) 


< C 6 £r 


(5.29) 


The assumptions (5.24) easily show that {pn) is uniformly bounded below by a positive 
constant, and it is then easy to deduce from (5.28) and (5.29) that the assumptions (5.17) 
of Theorem 5.4 are verified, provided we chose 5 small enough. Moreover (5.29) shows 
that £n and 7,1 differ only by a multiplicative constant that is bounded above and below. 
Consequently, in the notations of Theorem 5.4, w.h.p., we have 


.tn 


(n) 


m 


(n) 


-(5^ n7„ < |ei(n,^("^)| < 


m 


(n) 


m 


(n) 


+ 5 n'Jn, 


|e2(ra,C^"^)| < Sn-fn, 


and 


- 'Ipn 

njn 



Sn 


<< 5 , 

2 


-Wn < S. 
n 


This already proves the very last part of the result. On the other hand, using (5.28) and 
(5.29) allows to replace, up to a factor C5£n, by pnm\^^ and 7 ^ by £n, 

and a couple lines of computation provide the first part of the statement. To conclude, 
notice that (5.27) ensures that 


®n-©n 


< 5 

3 
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w.h.p. Hence, up to a factor 5, we may replace, in the above equation, by 


= PnXyG'^{{pnX + I - pn)y) 


—)^y^n + 1- — + en—^{x - 1) y + 0(en) 

mX \ \^2 ”^2 '^2 / / 


TTi 

= -^^yGn{Qn{x,y)) + 0{en), 


(n) 


m. 


(n) 


where both lines are easy Taylor expansions and O(-) is uniform in the parameters. This 
readily allows to conclude. □ 

We finally just need to prove Lemma 5.7. 

Proof of Lemma 5. 7. Let us prove the first inequality, the second being similar. First split 
the sum on r < l/on and r > l/on- To deal with the second sum, just write 


k=0 


k=l 


to obtain 


Ef 


r>llan k=0 


< ^2 r'^ 
r>llan 

<2al+^ Y, r^+VW(r) 

r>l/a„ 

< 2Mal+2 


so that Markov’s inequality yields 


\r>l/an k=0 


AM 


1+5 


(5.30) 


For the r < l/a„ term, first note that, for each r, the vector {nf,p^\k,r))o<k<r is 
multinomial with parameters np^^\r) and {[p]k)o<k<r- In particular, nf,p^\k,r) has a 
binomial distribution with parameters np^'^\r) and [p]^, so that Chebychev’s inequality 
yields 


Ec 


nilp-Kk^r) - np^^\r)[p\l 1 < y^n/iW(r)[p]^(l - \p]l) < A/n/rW(r)[p]^. 
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We can then compute, for a constant C depending only on M and r], that 



r 

E E^^ 





r<llan k=0 


r<llan k=0 



r<l/an k=0 




r<l/an 


k=0 k=0 


< C-^ \l(r) 

Jn ^ * 


r<l/an 


< c- 


y/nar. 


r2y/iW(r) 


r<l/an 


1/2 


< c- 




^/ho; 


\r<l/an 


E ^ 

\r<l/an 


-{i+v) 


< C 




where we use Cauchy-Schwarz twice. We conclude with Markov’s inequality that 


\r<l/an k=0 


> 2 “^ I - 


2C 1 


S yjna^ 


what, along with (5.30), provides the result. The proof of the second statement is similar, 
by merely writing in the second step of the reasoning. □ 


6 Gelation in a dynamic configuration model 

6.1 Introduction 

As already explained in Section 2.3, we can compare our model to a dynamic CM, at 
least until the first gelation time. More generally, we will introduce in the next section an 
alternative model that will allow us to compare our model to a DCM at any time, at least 
for the particles in solution. In particular, we will have to consider an arbitrary number 
n of particles with a degree distribution /i, where n < N and /r can be any (reasonable) 
distribution. However, the threshold for the size of the components that fall into the gel 
remains unchanged, still a{N). 

So let us consider a DCM described as follows. Consider n vertices, whose degree 
distribution is given by /r, and let 9 = (n, fi). Set i.i.d. clocks on the arms with distribution 
£(1), and link the first arm that rings with the second one, the third with the fourth, and 
so on. This therefore provides a graph process {9t{9)). If, at some time, there is only one 
activated arm, we assume that it does not matter in the graph structure. As usual, we 
define mf as the factorial moments of //, as its generating function, and we introduce 

= fex + 1 - y, F^{y) = yG'^{y). 
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As in the previous section, we will take constants c,m,M £ (0,+oo), and we will 
assume that 6 £ which was defined in Section 4. We will then obtain results with 

a precision k > 0. The results will only depend on these quantities, and we classically 
write K = K(c, m, M, k) to mean that A is a constant that depends only on c, m, M and 
K. Our main assumption on (a{N)) will be Assumption 2.3, and this will remain fixed, 
so we will not include it as a parameter. We will write 


“ a(iVp72 + N ) 


( 6 . 1 ) 


and all our results will be obtained with probability at least 1 — Kpi^, what we will shorten 
as w.h.p. in the proofs. 

Because 9 £ the generating functions have a fifth moment that is uniformly 

bounded above. In particular, if we do a Taylor expansion of (up to the fourth order), 
the constants hidden in the O(-) can be chosen uniformly in all the fixed parameters, and 
similarly for expansions of up to the third order, and so on. Similarly, the bounds on 
the will readily ensure that all the expressions involving the that we will consider 
are uniformly bounded above. We will implicitly always consider N large enough such 
that everything that we write makes sense. Specifically, since we assume (2.2), then any 
expression that we will consider multiplied by a{N)/N will be small, for large enough N. 
For the same reason, this N can be chosen uniform, and we can change the probabilities in 
the proofs to make the statements trivial for these small N. We are then able to perform 
Taylor expansions at the first order in a{N)/N (or a{N)/n), and whatever constant 
hidden inside the O(-) will depend only on the fixed parameters c,m,M. 


6.2 Evolution of the random graph 


Before diving into the heart of our problem, let us study the activated arms in the graph 
9t{9). At each time t, each arm is activated independently with probability pt = 1 — e“*. 
The law of large numbers then readily tells approximately how many arms are activated 
on the vertices of degree r, how many activated arms there are in total, and so on. To be 
precise, we will be interested in the following quantities. 


• ^k,r{9,t) is the number of {k,r) vertices at time t; 

• <5o{t) is its generating function; 

• ^-(1,1) is the number of activated arms at time t; 

ox 

• ?(0, B) is the first time t when ^e{t) = B. 

Finally, we let qt{x, y) = {ptx + 1 - pt)y. 


Lemma 6.1. Assume that (2.2) holds, and consider 6 = (n,/i) £ m and 5 > 0. 
Then, for some constant K = K{c,m, M, k), with probability at least 1 — Kp^, it holds 
that 

and 

-*Be(t) -ptmf 
n 


< 5 


a{N) 
N ' 


< 6 


a{N) 

N 
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In particular, for any B G [0, (1 — m)min], 


?(6»,S) +log (l -^ 

V nmf 


< 5 


a{N) 
N ' 


Note that the last quantity is only relevant when B £ [0, nrui], and we need to assume 
furthermore that B/{nm^) is bounded away from 1 to allow uniform Taylor expansions. 

Proof. This is again Lemma 5.7 with an = a{N)/N, and p = pt- It implies that 


r>l k=0 




< 5 


a{N) 

N 


w.h.p. The first part is then a simple computation and the second and third part follow 
easily. □ 


6.3 Gelation time and properties of the gel 

The following result concerns the time of the appearance of the large cluster, and its 
structure. It will be our main tool to prove precise estimates on the gelation times in our 
model, and on the evolution of the characteristics of the particles remaining in solution. 
We consider the following quantities. 

• is the time when a large component, i.e. a component of size > a{N) 
appears; 

• s^^\6) is the size of that large component; 

• (9) is the number of {k,r) vertices in that large component, and Us gener¬ 
ating function; 

• •) is the generating function of the number of vertices with degree r 
in that large component; 


6W(0) = 




(N) 


dx 


-(1,1) is the number of activated arms in that large component; 


Proposition 6.2. Assume that (2.2) holds, and consider 9 = (n,/i) £ m ^ > 
0. Then, for some constant K = K(c,m, M, k), with probability at least 1 — Kp^, the 
following statements hold. 

1. The gelation time verifies 


uW(0) + log 1-^ - 




run 


a{N)^ __ 

n 2 m 2 (m^ — mf) 


< K 


a{N) 

N 


2. In particular, the eharacteristics of the partieles at this time is deseribed by 




< K 


a{N) 
N ' 


and the total number of aetivated arms verifies 




n 


rrin 


n 2 {m 2 Y 


< K 


a{N) 
N ' 
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3. The size of the large eomponent enjoys 


a{N) < s^^\e) < {l + K)a{N). 


4- The strueture of the partieles in the large eomponent is deseribed by 


In partieular, 


and 


5. Finally 


1 ,(7V) 1 


a{N) 




m 


_ J^F 

^ mf ^ 


< K. 


a{N) 
1 


< K 


b^^\e)-2 


a{N) 


< K. 


< K. 


Proof. Let 

In the following, we shall consider of the form 


- m 1 + M 1 

= ^- sup a{k). 
4m ^ n n<k<nlc n n<k<nlc 


_ a{N) mf 

— ^n,± — o M M V-*- ^ 

n 


where we can assume k e (0,1/2). According to the moment assumptions, e~ < En < E'f. 
Moreover, since {a{N)) is such that (2.2) holds, then (ej/) verifies Assumptions (5.23). As 
usual, C will be a constant that might change from line to line but only depends on the 
fixed parameters, and we let 5 = njiflC) and 


af\«) = - log ( 1 - ) - log (1 - (1 ± J) 


m. 


2m^( 


mo — va\ 


We may then take K as in Theorem 5.6, which depends only on d, r], m, M and (e^), 
and thus only on k, r], m, M, c and (a(A^)). All the conclusions of the corollary will thus 
hold with probability at least 


1-K 



and thus w.h.p. 

1. Denote 

Pn =P^W^e) = ^{^+en,±). 

Now, note that at any time, St(0) is exactly a CM{n, graph. Therefore, > 

a^\6) means that CM{n, y,Pn) has no large component. But Theorem 5.6 implies that. 
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w.h.p. 


ei(n,/i,p+) 


\ m^mg J n 


= a{N) 1^1 + K — (5(1 + k) 

> a{N) (1 + K - C6) 

= a(iV) (i + k/2). 


m 


M 


2m I 




Hence, < a^\6) w.h.p., and similarly, a^^\9) > a^\9) w.h.p. Then, an easy 


(N). 


Taylor expansion allows to write 


(t^±\9) = - log ( 1 - ^ ) + (1 ± (5) 


mi 


a{N) 


mn 


mn 


L2 ! n 2m2{m2 — m^) 

and Point 1 of the result follows. 

2. Since is increasing in t, from the point above, w.h.p. 

<&e{^^-\9)) < <5g{<j^^\9)) < ©0(af V)) 

On the other hand, by Lemma 6.1, w.h.p. 


+ 0 


»{N) Y 
N J ’ 


\9)) < + 

But a simple Taylor expansion shows that 

Another Taylor expansion ensures that 

G'/.(g^ov)(g^(x,y)) = Gf,{qf,{x,y)) + “ '^)yG'^i<l^^{x, y)) + ko(^^^ 


N 


But these computations actualy hold in the || • ||2 sense, as well as for (5g{a[!^\9)), and 
the result for <3g{a^^\9)) follows. The second part is then a simple computation. 

3. Note that the size of the largest component in 9t{9) is increasing with t. But we just 
saw that w.h.p., a^^\9) < a^\9), so that 

a^^\9) < ei(n,Ai,p;! 

w.h.p. A similar computation as in Point 1 above ensures that this size is bounded by 
1 + K + G5 < 1 + 2k w.h.p. 

4. Denote by ^jJ^ the generating function of (ufc^r(Ci(n, fj,,p^))). Theorem 5.6 implies that 
w.h.p.. 


1 


nSr, 


/ / to 


< 6 . 


Rearranging, it is easy to see that 

1 


a{N) 


'^B - ^xyG'^iq^) 

11 L-\ 


< Gd < K. 
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A similar inequality holds for . But 

Vk,r{Gi{n,lJ.,pp)) < <Vk,r (ei(n,/i,p+)) 

w.h.p., so that 




a{N) 


m‘ 


< n. 


The other formulas are obtained after easy computations. 

5. Finally, for the same reason, 

t'u;W(6») <- ^^^fc,r(ei(n,/i,p+)) < K 

Ti ‘ ^ 


n 


k,r>0 


W.h.p., and the proof is over. 


□ 


6.4 After gelation 


The previous result precisely describes when gelation occurs and what happens for the gel 
at this time. It therefore allows us to explain what remains in solution after a gelation 
event. Let us then start by defining the following quantities. Note that the notation is 
the same as above, but with capital letters. 


• S^^\8) = n — s^^\9) is the number of vertices remaining in solution at 

• — ^k,r{8,a^^\6)) — is the number of {k,r) vertices in solution at 

and is its generating function; 

• •) is the generating function of the number of vertices with degree r 
remaining in solution at 


• B^^\9) 
a^^\9). 




(1,1) is the number of activated arms remaining in solution at 


We then have the following result, which is a direct consequence of Lemma 6.1 and Propo¬ 
sition 6.2. 

Lemma 6.3. Assume that (2.2) holds, and consider 9 = {n,pL) G M k > 0. Then, 
for some constant K = K{c,m, M, k), with probability at least 1 — Kp]\f, the following 
statements hold. 


1. The number of particle remaining in solution enjoys 

n - (1 + K)a(iV) < S^^\9) <n- a{N). 


2. The structure of the particles remaining in solution is described by 

if'")(«) - GM - - T.vG'„m 


< K 


a{N) 
N ' 


3. In particular. 


and 


n 


a{N) 


+ 

n 

iiFu ^ 

2 

N 


a{N) / 


< K 

a{N) 

n 1 

2(m^)2 y 

N 
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6.5 Properties of the new model 


After gelation, we are left with particles in solution described by the above result. In par¬ 
ticular, there are S^^\9) of them, and their total degree is described by It thus 

makes sense to ask what happens for a configuration model with these parameters. It will 
turn out to be especially useful when we consider the alternative model in the next section. 
To be precise, we consider S^^\9) particles, with degree distribution ]1^^\9), which has 
generating function S^^\9). We write 9^^'^ = {S^^\9),]J^^\9)). The 

factorial moments of JI^^\9) are written m\^\9), and we denote the time when a giant 
component appears in as a^^\9). We can then prove the following result. 

Lemma 6.4. Assume that (2.2) holds, and consider 9 = {n,^) G m k > 0. Then, 
for some constant K = K{c,m, M, k), with probability at least 1 — Kp^, the following 
statements hold. 

1. The generating function of Jl^^\9) verifies 




a{N) 


n 


C _ —F 

fiL /i M 


< K 


a{N) 


N 


2. The factorial moments of pL^^\9) enjoy 

m] [9) — mf - mf - u — 

n \ mf I 


and 


and 


’{9) -mif- I 


< K 


< K 


a{N) 
N ’ 

a(A^) 
N ’ 


m^^^\9) — rug 


< K. 


7)W^ 


3. The time when a giant component appears in S{9^ ) verifies 


aW(0) + log(l-^)- 


a{N) 


Srug 


1 


m 


2 I n I 2772-2 — m)') 


< K 


a{N) 
N ' 


Proof. 1. We know from the previous result that w.h.p. 


(TV) ^ a>W(0) 

^ Si^){9) 

1 1 
“ SW(0)/nn 


$(^)(0) 


1 


1 — a{N)/n + K 0{a{N)/n) 


Ga- 


a{N) 1 


Fa ~\- K O 


n m‘ 


1 


= 1 + 


a{N) 


n 


+ KO{a{N)/n)^ (^G^- 


a{N) 1 


LL^ fd, 

n nil 


Fa F K O 


a{N) 

n 

a{N) 

n 


— Git + 


a{N) 


n 


Ga - JlFa ] + K O 


m 


1 


a{N) 


n 


and this holds in the || • ||2 sense. The first part of the result follows. The estimations of 
m[^\9) = {Gg^^y{l) and m^^\9) = (Gg^^)"(l) are then easy computations. 
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2. For the third moment, note that by Proposition 6.2, 


= 




sW{e) 


n 


SW{6) 




3.,W/ 


fc,r>0 


fc,r>0 


m 




n 




= (1 + K 0(1)) (m|( + K 0(1)) 
= TTig + K 0(1), 


and the result follows. 

3. Finally, the estimation of S^^\6) and m\^\9) above readily shows that the assump¬ 
tions of Proposition 6.2 are still in force, maybe at the price of considering different 
constants. Therefore, Proposition 6.2 applies, and shows that w.h.p. 


7fW(a\ = - 




a 


(0) = - log 1 - 


m 


1 


ie) 




+ 


a{N) 




m 


io) 


2 {9) J 2m2^\9){m2^\9) — m[^\9)) 

A couple lines of computation using the previous results show that 


+ kO 


a{N) 

N 


m[^\9) 

rn^2^\9) 



m2 

{rriiY 


+ 




m'^rrio 


+ K 0 ( 1 ) 


A couple more lines ensure that 


log 


m[^\9)\ 
m2^\9)) 


log 1 


K I ^ «(^) I ^ 

1712 j ^ \mi 


m 

1^2 {mi^ 


3 



+ kO 



The corrections on the second term of a^^\9) are of the second order, and thus, one may 
remove the bars there at the price of another factor k 0{a{N)/n). An easy computation 
then allows to conclude. 

□ 


6.6 Non-degenerate conditioning 

Recall from Lemma 4.7 that, right after a gelation event, conditionally on the particles 
remaining in solution and the number of activated arms, the structure in solution is that 
of a CM on these particles with this number of activated arms, conditioned on having no 
large component. It therefore makes sense to ask whether this conditioning matters or 
not. The main interest of the previous result is that, w.h.p., this conditioning does not 
matter. It will be an easy consequence of the following computation, which we will reuse 
later. 

Lemma 6.5. Assume that (2.2) holds, and consider 9 = (n,/r) E M k > 0. 
Then, for some constant K = K{c, m, M, k), with probability at least 1 — Kp]\[, 




a{N) / rrig 1 \ ^ a{N) 

n y2m2 (m ,2 — ruj') ruj* y “ N 
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In particular, 




a{N) 

mf 

n 

7772 ( 77^2 — m^) 


< K 


a{N) 
N ' 


Proof. From Lemma 6.4, we know S^^\6) and m\^\9) with good precision. 

Easy computations allow to obtain 


B^^\9) _ K A a(A^) 


rrin 


n 


2m f 


mo 


+ 


— m^ 

{mfY 


+ kO 


a{N) 

N 


Therefore, w.h.p., and maybe at the price of taking smaller constants, 

B^^\9)/{S^^\9)m[^\9)) < 1 - m, 

and, as we mentioned in the previous proof, 9 G We may then use Lemma 6.1 to 

write 


c,{9,B^^\9)) = - log (^l- 




S(^){9)m[^\9), 


+ kO 


a{N) 

N 


and the result follows after easy computations using Lemma 6.3. The second part is then 
a direct consequence of Point 3 of Lemma 6.3. □ 

Proposition 6.6. Assume that (2.2) holds, and consider 9 = (n, /i) G m ^ > 
0. Then, for some constant K = K(c,m, M, n), conditionally on 9^^^ and B^^\9), a 
configuration model CM{9^^\ B^^\9)) has no large eomponent with probability at least 
1 - iLpAT. 

Proof. From the previous result, w.h.p. 

rW(0) _ ^(0^ B^^\9)) > 






n \ 771,2 (m ,2 “ 


— K 


The fact that 9 G m ensures that, if k is chosen small enough, then this quantity 


is positive. Therefore, w.h.p., gelation in a occurs strictly after the time when 


B(-^\9) 


arms are 


activated, and thus a CM{9^^\ B^^\9)) has w.h.p. no large component. 

□ 


7 Alternative model 

After these two sections concerning random graphs, we come back to our model, and start 
putting the pieces together. The goal of this section is to present an alternative model 
which is easier to study than ours since it has a much nicer combinatorial structure. We 
will prove however that w.h.p. the two models have the same distribution on compact 
time intervals. 

According to Lemma 4.8, at any gelation time r and conditionally on the number 
of particles in solution n = |5'(t)|, their empirical measure of degrees fi = fir and their 
number of activated arms B = B(t), the configuration in solution is that of conditioned 
configuration model CM'{n, fi, B). Informally speaking, our goal is to somehow get rid of 
the conditioning. 
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We introduce below a slightly more complicated process, which is however merely built 
from DCM. Encoded in that process are two processes: our original one (at least, its state 
in solution), and another one which can be more easily studied. In other words, we have 
a coupling between our process and a more simple one: this is the purely combinatorial 
Lemma 7.2. The fundamental fact is that the processes are actually equal w.h.p., as 
proved in Lemma 7.3. This second statement makes heavy use of what was proven in 
Proposition 6.6. This construction and the results are very similar to those presented in 
Section 5 of [23] , to which we refer for a more elementary exposition. 

For each > 1, consider a family k)) of independent DCM, indexed hy k >1 

and 0 = with n G [A^j and fi G (these were dehned in Section 4). Note that 

there is only a countable number of parameters and there is thus no issue considering this 
family. For each such 6 and k, %^^\9,k) has the distribution of S(0). We shah use the 
same notation as in Section 5, but adding a parameter k to insist that we are dealing with 
S^^\0,k). Specihcahy, we let 

• R, k) to be the hrst time when there are R activated arms in k)] 

• a^^\9, k) be the time when a component of size > a{N) appears in k); 

• s^^\9, k) the size of that large component; 

• B^^\9, k) the number of activated arms remaining in solution at a^^\9, k); 

• ]I^^\9, k) the empirical distribution of degree of the particles remaining on solution; 

• 9^^\k) = {S^^\9,k),-p^^\9,k)). 

Similarly as in [23], let us construct the process {Y{t)) as follows. 

Step 0 Let po = to = 0, Dq = 0, and K(0) = 1. Consider 2^^\pq,K{0)), 

up to the time ti = a{po, K{0)) when a large component appears. Then let pi = 
p^(^)(K{0)), Di = B^^\po,K{0)), dehne y(s) = si^^(po,^^(0)) for s G [0,ti), and 
go to step 1. 

Step i Consider here the graph processes 9^^\pi,k) for A: > 1. Dehne K{i) to be the 
hrst k such that n k) has no large component. Then take 


ti+i =ti + a^^\pi, K{i)) - Di, K{i)), 

as well as pi+i = p^^^\K{i)), A+i = B^^\pi,K{i)), dehne 


F(s) 


qW 

(pi,Di,k)-ti-lrS 


{pi,K{i)), 


s G [tj, 


(7.1) 


(7.2) 


and go to step i + 1. 

Obviously, we stop and let tj+i = +oo when no more large component can be created. 

An important thing to notice, and the main difference with [23], are Equations (7.1) 
and (7.2). To understand where it comes from, let us explain the construction in more 
details. At Step 0, we are just looking at a DCM 9^^\po, 1) until the gelation time ti, just 
like our original process. Thanks to Lemma 4.8, we know that, right after gelation, what 
remains in solution has the distribution of a CM'{pi, Di). Now, we want to replicate this 
latter graph, thanks to some DCM k). To this end, we look at these processes up 

to the time (;^^\pi,Di,k) when Di arms are activated. At this time, we obtain graphs 
with distribution CM{pi,Di). We choose the Ar(l)-th, the first one which has no giant 
component, and which has therefore the conditioned distribution CM'{pi,Di). By the 
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Markov property, the evolution of our original process and of afterwards 

are similar. 

In particular, note that we are dealing with two different time lines here: the one of 
and the one of K{1)). We look at the first process on the interval 

[0,ti), and at the second one on the interval [c^^^\pi,Di,k),a{pi,K{i))). Therefore, 
Equations (7.1) and (7.2) merely allow to define our new process Y by pasting these two 
parts one right after the other. In the definition of [23], the times were automatically 
matching and no such treatment was necessary. 

Remark 7.1. It might seem more natural to eondition on Nt and tti and use Lemma 
4-9. However, the whole triek of the previous eonstruetion is to he able to replieate a 
CM'{pi,Di) thanks to a DCM, so as to he able to study the dynamie afterwards. As we 
saw, it suffiees to stop at the right time. But, even replieating a CM{Nt,'Kt) does not seem 
possible, sinee there is possibly no time when we get exaetly the distribution of aetivated 
arms vr^. This explains the slightly bizarre eonditioning on the total number of aetivated 
arms. 

By iterating this reasoning, it should be clear that we have constructed a process which 
evolves in the same way as our original process in solution, at least in between two gelation 
times. Our construction does not ensure continuity (in terms of graph structure) at the 
gelation events, but it does not matter for all intents and purposes. Doing so is possible, 
but would further complicate the definition. Therefore, with a slight abuse of notation, 
for two processes U, V and two sequences of random times (oj), (6j), we define 

if 

((ai),C/(a,)))^=^( 60 ,^( 6 d))> 

where, for an interval I, Uj = {U{t),t G I) and ^ 0). In the following 

result, recall that {Gs{t)) is the configuration in solution of our process and (ri)i>i is the 
sequence of gelation times. We also let tq = 0. 

Lemma 7.2. The equality 

((n),Gs)‘='{{t.),F) 

holds. 


Proof. Let us write the parameters of our original model as Bi = Bt-^ 

and Ci = (G5)[ri,Ti+i)- particular = Po- Now, our process is, before 

gelation, just a S(^o)) so has the distribution of 9^^\po, !)• In particular, the time of the 
appearance of a large component and the properties of what remains in solution are the 
same as for Y, i.e. 

Let us then prove by induction that for every z > 1, Ddj holds, where IKj is the assumption 


l<j<i) = {[tj-i,tj,pj,Dj,Y[t._^^t.^],l <j< i). 


We just checked “Ki, so assume that “Ki holds for some z > 1. On the one hand, we know 
from Lemma 4.8 that 


Gs{Ti) GM'iCuBi). 
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Now, from the construction, V(ti) is constructed from the graphs ^)- Each 

has distribution CM{pi, Di), and we choose the K(i)-ih graph, the first which has no 
large component. It is thus distributed as a CM'{pi, Di) graph. By dfj, this has the same 
distribution as a CM'(^i, Bi) graph, and thus 

Y{ti) CM'{pi,Di) Gs{n). 

Then, by Markov property, Gsiu + t)t>o and Y{ti + t)t>o evolve in the same way until 
the next gelation event, and TCj+i follows. □ 

Hence, from now on, we can assume that (Gs{t)) is constructed as above. In particular, 
it is coupled with the process that we introduce now. We write !)• 

Step 0 Let 6q = co = 0 and Aq = 0. Consider up to the time ui = 

a{0Q, 1) when a large component appears. Then let 6i = 1), 

define Z{t) = 9i^\do) for t G [0,(Ti), and go to step 1. 

Step i Consider here the graph processes Then take 

C7i+i = a* + VO, (7.3) 

as well as = 6i^^\l), Hj+i = B^^\9i, 1), and define 

Z{s) = 1)> ^ e ko (7.4) 

where this interval might be empty if fii+i = cjj. Then go to step i + 1. 

In other words, we consider the same model as above, but taking always K{i) = 1. It 
should be clear that {Z{t)) is easier to study than (Y(t)). Our next claim is that under 
(2.2), Z and Y are barely any different, at least on compact intervals. To this end, define 

I = inf{i > 1, K{i) + 1}, E{T) = {ti > T}, 

so that E{T) is the event that K{i) = 1 for all the 9^^^ we consider before time T. Hence, 
on E{T), the processes {Y{t),t G [0,T]) and {Z{t),t G [0,T]) are equal. As usual, we 
will write E^^\T) when we want to insist on the dependence on N. We will prove the 
following. 

Lemma 7.3. If Assumptions 2.2 and 2.3 hold, then, for all T >0, F{E^^\T)) —^ 1. 


Proof. Take c, m, M as in Lemma 4.1. We may assume that m is small enough such that 


- log 1 - 


1 


1 + m 


a{N) 

N 


> T. 


(7.5) 


Consider k = 1 for instance, and the corresponding K = K(c,m/2, M, n) such that all the 
results of Section 6 hold. Note that we take the constant corresponding to m/2 to leave 
ourselves a bit of leeway, what we will use to study F^^\T,2) below. All that we write 
from now on is implied to happen with probability at least 1 — Kp^, where we recall that 
Pn is defined in Equation (6.1). 

We want to prove that F{E^^'>{T)) —>■ 1. Define 

J = inf{i >l,pi^ 

and the event 

EW(r,l) = {vt G [0,r] G 

considered in Lemma 4.1. There are three ways things can go wrong, i.e. an outcome oj 
is such that tj{uj) < T: 
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• either does not occur; 

• or E^^\T, 1) occurs, but J < I and tj <T; 

• or E^^\T, 1) occurs, J > I, but K{i) ^ 1 for some i < I. 

Let us call these events respectively 1) = E^^\T, l)'^, E^^\T,2) and E^^\T,3). 

By choice of the constants c, m, M as in Lemma 4.1, we already know that 

F{E^^\T,1)) 0 . 


Let us now consider E^^\T,2). In this case, since tj < tj < T and E^^\T,1) 
occurs, then pj = {nj,pj) is in ™ i.e. it verifies nj > cN, 

SrM < M, mg-^ > m, but m^'^/rrig'^ < 1 + m. Moreover pj-i £ 

But recall Lemma 6.4. It easily implies that there is a constant C such that, for any 
p = {n,p) eK[^2,M^ 

mfV) ■ ^2 I N )• 

Since m^/mg > 1 + m, then m^'^^(p)/m 2 ^^(p) > 1 + m/2 for N large enough. Therefore, 
since pj-i G then 






> 1 + m/2. 


Hence pj G 
gelation time verifies 


According to Proposition 6.2 (with k = 1), the corresponding 


>-log h - - 


a{N) 

N 


As we mentioned, this is our only reason to choose the constant K corresponding to m/2, 
not m. Now, since < 1 -|- m as we just mentioned, then a^^\pj) > T by the 

choice of m as in (7.5), so that tj > T. To conclude, P(F(^)(T, 2)) < 1 — Kpj^ —)■ 0. 

Finally, on E^^\T, 3), all the pi we consider are in Note that K{i) = 1 means 

that 9^22 D large component. But 


SS,D.,Ai)(P-l) = CM{pi,Di) = CM(^(^),hW(p,_i)) 

and we know by Proposition 6.6 that this graph has no large component with probability 
greater than 1 — Apjv- Since there are at most N/a{N) gelation events, and recalling the 
definition of pN, we get 


P(F<»l(T,3))</r^P. 





0 


by (2.2). The result follows. 


□ 


Notice that in the very last computation we crucially use (2.2). This is the main reason 
for our assumption, just as in [23], we needed to assume slightly more on the threshold 
{a{N)). As a by-product of the above proof, namely the fact that P(F('^)(T, 2)) —)■ 0, we 
see that all the (n*, pt) that we consider before time T are well-behaved, in that they are 
in k22,m- Coupled with Lemma 7.2, this allows to obtain the equivalence of our model 
and the alternative model on any compact interval, along with the good behavior of the 
relevant quantities. 
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Theorem 7.4. If Assumptions 2.2 and 2.3 hold, then, for allT > 0 
¥{{Gs{t),t€ [o,r]) = {z{t),te [o,r])) ^ i. 


Moreover, there exists c,m,M G (0, +oo) sueh that 


8 Limit 

This section puts together the results on the DCM and the alternative model to study 
precisely how our model evolves between gelation events. This culminates with the proof 
of the convergence of and the explicit characterization of its limit. 


8.1 Microscopic evolution 

Let us recall that Pt{k, r) is the number of particles in solution at time t, with k activated 
arms and degree r, and is the empirical distribution of degrees in solution at time t. 
Define 

1 _ ^ AT 

Mx,y) = ^ 

r>0k=0 

Thanks to the alternative model and the results of Section 6, we are able to describe 
precisely how ■i/’t and (ft evolve at the microscopic scale. 

Recall that (rfc) is the sequence of gelation times, and with a slight abuse of notation, 
denote p,k = hr,,, 'rnf = mf*" and (fk = (frk- Note that all these quantities are 

constant on [Tfc,rA;+i). We finally let V'fe = V’r” > ^he left limit of xf at Tk+i. 

In all this section, we will use freely the alternative model and Theorem 7.4. The most 
important point to remember is that, with high probability, our original model and the 
alternative model coincide before time T, all the 9k we consider before this time are in 
^2rn and all the K{i) are equal to 1. In particular, our model on the interval [r^, Tk+i) 
corresponds to a DCM with parameter 9k on an interval [‘;^^\9k,Bk,l),(T^^\9k,l)). In 
all the proofs, O(-) is implied to be independent the k that we consider, and we imply 
that everything happens w.h.p. 


Lemma 8.1. Assume that Assumptions 2.2 and 2.3 hold. Then, the following hold. For 
any k > 0 and T > 0, there is a K > 0 sueh that, with high probability, for all k > 1 sueh 
that Tk < T, we have 


a{N) 


m 


k-l 

3 


m 


k-l 

2 


(m 


fc-i 

2 


— m^ 




) 



< Tk+l 


Tk < 


a{N) 


m 


k-l 




k-l 


k-l\ 



and 


and 


(fk+l 


(fk + 


a{N) 


1 


m 


k-l 

1 


y4>'kiy) 


< K 

2 


a{N) 

N 



4>k 



X + 1 
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In particular, for any T > 0, there are constants d,D ^ (0, +oo) such that, with high 
probability. 


n(N) 


Tk < D 


a{N) 

N 


for all k > 1 such that Tk < T. 


Proof. A first consequence of the alternative model is that we have 


r.+i - r, = B,) = (4'!>. . 

and from Lemma 6.5, this is 

aW 1 "4'‘ 

Nt V ~ L 

The second consequence is that 


f’k+l 


N 


and Lemma 6.4 shows that this is 

Nk+i a{N) 

N \ Nk 


nil 


By Lemma 6.3, A^fc+i = — a{N) + kO( 1), and it is then an easy computation to 

conclude. 

As a third consequence, we obtain 



According to Proposition 6.2, and this is 


G^ki^Uk) + O 




X + l 




Since there are less than N/a{N) gelation events, these relations are true for all k with 
Tk < T, with probability at least 1 — KNp]sf/a{N), and as we saw, this tends to one 
because of (2.2). 

The last part of the result is a direct consequence of the first inequality and the fact 
that all the 9^ that we consider are in w.h.p. □ 


8.2 Tightness and limits 

The result above allows to describe precisely the evolution of (r^), {(fk) and (V’fc)- It 
is therefore not surprising that it allows to describe limits. Classically, we first prove 
tightness, then show that subsequential limits solve an equation, and finally solve this 
equation. Recall that the space is the space of k times continuously differentiable 
functions on [0,1] or [0,1]^ (the context should make obvious which one) endowed with 
the II • llfc norm defined in (2.10). 

Lemma 8.2. Assume that Assumptions 2.2 and 2.3 hold. Then the sequence is 

tight in D(M'’',M), and any limiting point is locally Lipschitz-continuous. Moreover, the 
sequence is tight in D(]R+, C®), and any limit point is continuous. 
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Proof. Concerning and this can be done exactly as in [23], using the classical 

criterion of [7]. From the previous result, the number of gelation events on an interval of 
size s is at most a{N)/Ns/d, and each makes (resp. in || • jjo norm) change by 
at most 2a{N)/N, see [23] for details. 

Note that this only shows tightness of in However, this implies the 

tightness of in so in ^'^(N) thanks to Lemma 4.4, and this, along with 

the tightness of implies in turn tightness of in D(M+, C^). □ 

Proposition 8.3. Assume that Assumptions 2.2 and 2.3 hold. Then any limit point (n, (p) 
of mD(]R.+ ,]R.) X D(]R.“*“, C®) verifies 


— 1) (pt — ^ — ^gel) 


and 


where 


4 ,., == 


y(p't{v), t > Tgel, 


s{t) 


1 cP'l{l){cP'l{l) - cP'm 

nt 0'/(l) 


Moreover, whenever eonverges on a subsequence m C®), so does {^p^^'^), and 

its limit Ip then verifies 


Mx,y) 



X + 1 



Proof. The part of the result before Tgei is clear since we know that w.h.p., no gelation 
event has occurred before Tgei (thanks e.g. to (2.6)), and therefore nt and <pt are constant 
on [0,Tgei). We extend to Tgei by continuity. The formula for (pt is then a consequence of 
Lemma 6.1. 

For the part after Tgei, hx Tgei < t < t + s, and define 


r^I^\t,t + s) = 


= inf 


m 


i-^u 





{t,t + s) 



All the following events will happen w.h.p. First, at least one gelation event has happened 
by time t by (2.6). Moreover, by Lemma 8.1 and by monotonicity of {Nt), we have, if 
t < Ti < t + s, 


a{N) 


(r^.^\t,t + s) - 


< Tj+l -Ti < 


a(N) 


Nt 


t+S 


(r^_^\t,t + s) + HiJ . 


The number of gelation events on [t, t + s] is thus at least 


Nt+s 1 

a{N) r^^(t, t + s) + K 


and at most 

Nt 1 

a{N) r^'^{t,t + s) — K 

By Proposition 6.2, a gelation event makes at most (1 + K)a{N) particles falls into the 
gel. Consequently 


Nt+s 1 

^ r^_^\t, t + s) + K 


L ''-t+s 


— n- 


(N) 


t ^ 


< - 


a{N) Nt 1 

N N g'j _ f. 


(1 + k). 
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But ^ in D(M+,C75), and (/> is continuous (in time), so there is actually uniform 
converge on the compact sets (see [7]). Therefore 

t + s) —)• inf l/s{u), r^^(f, f + s) —)• sup l/s(rt) 


and, passing to the limit in the equality above then shows that 

1 1 


-ni+< 


sup„g[t_t+^] l/s{u) + K 


< nt+s -nt < -nt- 


inf„g[t^i+^] l/s(t6) - K 


(1 + k). 


Having k —?■ 0, then s —?■ 0 and using continuity allows to conclude. The computation for 
(j) is similar. 

Concerning assume that converges along some subsequence to some (j). 

For t > 0, take k^^\t) the k such that t G Clearly, ip is an increasing function 

in time on each interval ), and cannot change by more than a{N)/N (in || • ||o 

sense) at each gelation event. Therefore, for all t > 0, 


a{N) 

1pkW(t)-l - - i’kW{t)- 


By Lemma 8.1, the RHS is 


V'fcW(i)(x, y) = I I ^ ^ I ^ I + ^ ( 


a{N) 


N 




4w(i) 


( 1 ) 


x + l- ^ 




( 1 ) 


y \ +0 


a{N) 

N 


( 8 . 1 ) 


( 8 . 2 ) 


and a similar formula for the LHS. Again by uniform convergence on compact sets, we 
have (pt- On the other hand, —>■ 0 by Lemma 8.1, so 

that —t (pt for the same reason. Passing to the limit in (8.1) and (8.2) shows the 

result. This actually just shows pointwise convergence, but (8.1) and the lower bound on 
Tfc+i — Tfc of Lemma 8.1 easily ensure that the variations of tp are of order s on intervals 
of size s, what allows to conclude classically to tightness, as in Lemma 8.2. □ 


We can finally conclude to the convergence of all the relevant quantities by solving the 
previous equation. 


Proof of Theorem 2.13. First, assume that ut{x) is a differentiable function that verifies 

(iut{x) = —xq{t)u^{x), t > T, x G [0,1], (8-3) 

where g is a continuous function. This can be solved by method of characteristics. Indeed, 
define Q{t) = exp — f.^ q(s) ds. Then it is easy to check that ut{x/Q{t)) is constant, and 
therefore equal to urix), and thus (8.3) has a unique solution given by 


ut{x) = UT{xQ{t)). 

Now, let us define ut{x) = p[{x)/{xp'f{x)). On the one hand, we have 

f>'t{x)cPT{x) 


// N 1 4 >'t{x) 

u^{x) = - 

X Xpl{x) 


x‘^p'1 {xY 


(8.4) 


(8.5) 
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Let us define q{t) = —nts{t)/(p^{l). A simple computation shows that 


Q{t) 


ut{l) - 1 
u[{l)+ut{l) - 1‘ 


( 8 . 6 ) 


But on the other hand, using the PDE solved by 4> of Proposition 8.3 for T = Tgei, it is 
easy to check that u solves (8.3), so that ut = UT{xQ{t)). Plugging this in (8.6), we see 
that Q verifies 

a(f) = _9^ = _ UT{Q{t)) - 1 _ 

Q{t) Q{t)u'j,{Q{t)) + UT{Q{t)) - I 
This shows that R{t) = Q{t){l — UT{Q{t))) verifies 


R'{t) = t > Tgei. (8.7) 

But, by continuity, (I)t = Gfj,, so that 

^ xG"(x) xGUxY 

Therefore 

R{T) = Q{T){1 - UTiQim = -ut{1) = 

rrii 

which, with (8.7) shows that R(t) = for t > Tgei. By definition of R, this is just saying 
that Q{t) solves (2.14). 

Finally, since </> also solves (8.3), we get from (8.4) that 


4>t{x) = (f)T{xQ{t)) = G^{xQ{t)), 


( 8 . 8 ) 


whence we deduce the solution for nt = </>t(l) and V’t given in the previous result. □ 


9 Conclusion 


9.1 Typical clusters 


With all these tools in hand, proving Theorem 2.14 is an easy task. According to Theorem 
7.4, it suffices to do it for the alternative model. Recall that, at each time, the alternative 
model {Z{t)) is the realization of a DCM at some time. In part, for each t, Z{t) is a CM 
with some random parameters Mt and By Theorems 2.13 and 7.4, we have 


N 


nt, 


AN) 


TTt. 


The second convergence holds in Hence, (Aj'^^) verifies the assumptions of Proposition 
2.5 with limit tt*. Proposition 2.5 then applies for Z{t), conditionally on Mt and Xt^\ 
and it then suffices to de-condition. To be precise, one can use Skorokhod representation 
and prove the result as (2.16). The easy details are left to the reader. 
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9.2 Limiting concentrations 

Let us now explain Formula (1.8). We assume that we are in the supercritical case 
> rni- The reasoning is again done as in Section 2.2.2, and we consider the same 
notation. 

First, note that Q{t) lives in [0,1], and that any limit point c has to solve (1.7). This 
equation has a unique solution, and thus Q{t) —)• x as t —)• +oo. Hence, taking limits in 
(2.15), we see that irt — p, where p is a probability with generating function 


Gp{x) 


G^{cx) 

G^(c) ■ 


Similarly to the reasoning done in Section 2.2.2, we have 


c^^(0,m) = (v{e[^J{Qo)) = m) . 


According to Theorem 2.14, Ct^^(oo) should be a GWp^p. This is of course not rigorously 
true since we only have convergence on compacts sets. Notwithstanding, 

lim Q (v{e[yJ{oo)) = m) = Q (GWp,p = m) = — 3— -/3™“ V*™(m - 2), 

N^+cx) \ j't' j m — 1 


where, again, the last equality is an application of Dwass’ formula [13]. 

This finally explains the appearance of the term /3 in Formula (1.8), compared to 
Formula (1.6). In the latter case, there is no gelation, and typical clusters in the final 
state are subcritical Galton-Watson trees. In the former case, there is gelation, and typical 
clusters in the final state are critical Galton-Watson trees GWp^p. Their reproduction law 
is obtained from the original subcritical law p, by an exponential tilt, which turns it into 
a critical law. 


9.3 Typical particles in solution 

Theorems 2.13 also shows how typical particles in solution evolve. Gonsider a particle p 
uniformly chosen at some time t > Tgei. Then, in the asymptotics A" —>■ oo, 

• p is still in solution with probability nt = Gp{Q{t))] 

• p is in solution and has degree n with probability p{n)Q{t)'^ (the coefficient of p” in 

• knowing that p has degree n, it is still in solution at time t with probability Q{t)^. 

The last point suggests that a form of asymptotic independence holds throughout the 
whole process between the distinct arms of p. Everything indeed happens as if each of 
these arms causes, independently of the others, the fall of p into the gel by time t with 
probability l — Another way to see this is that the growth of the different subclusters 
attached to each arm of a given particle most likely do not interfere with one another, so 
that the particle falls into the gel when one (and only one) of this clusters reaches size 
about a{N). 
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9.4 Particular cases 

9.4.1 Poisson distribution 

Assume that /i is a Poisson distribution y(A) with parameter A > 1. Then one readily 
obtains 

Tgei = - log (l - • 

Moreover, = fJ-, and for t > Tgei, Q{t) = + 1/A, 

rit = exp — A + Ae“*^ , 

and TTt = p '.= 1P(1). Hence, in this case, the distribution in solution is essentially that of 
a configuration model on ritN particles with degree distribution p. A typical cluster in 
solution is a GWp, and has thus a stationary distribution. This should not be surprising, 
since the configuration model with a 1P(A) distribution corresponds to the Erdos-Renyi 
graph with parameter A. This agrees with the results of [23], where typical clusters in 
solution are such GWp trees post-gelation. 


9.4.2 Fixed degree 

Assume now p, = 5d for some d >3, that is, all particles initially have degree d. Then 

d-2 


Tgei = - log 


d-l 


Moreover, u = and for t > Tgei, Q{t) = e — l)/{d — 2) 


rit = 


d-l 
d-2^ 


and TTt = p := Bin(d, l/(d — 1)). Finally, p = Bin((i — l,l/(d — 1)). In this case, the 
distribution in solution is essentially that of a configuration model on ntN particles with 
degree distribution p. A typical cluster in solution is thus a GWp^p, and is again stationary. 

Notice that in this case, nt —)• 0, which makes sense since there are no particles with 
degree 1, so that, in some sense, clusters cannot “close” by attaching such particles at 
their extremity. 


9.4.3 Binomial distribution 


Assume finally that p = Bin((i,p) for d >3 and p > l/{d — 1). Then 


and IT = Bin((i — l,p). The computations are slightly more tedious in this case, but it is 
nonetheless easy to obtain that, for t > Tgei, 


nt = 


d-l 

d-2 


(1 


-p + pe *) 


d 


as well as TTt = p ■= Bin((i, l/{d — 1)) and p = Bin((i — 1, l/{d — 1)). 

We therefore have the same typical clusters in solution GWp^p as in the previous case. 
This is easy to explain. Note indeed that, before gelation, the number of activated of 
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the particles has distribution Bin((i, (1 — e“*)) in the previous case, and Bin(d,p(l — e“*)) 
in the current one. At gelation, typical clusters are therefore in both cases, and 

then the evolutions are similar. There is however a non-trivial time-change due to the 
difference in degrees. In the binomial case, there is moreover a positive proportion of 
particles with one arm, and thus there remains a positive mass at the end. In the final 
state, we essentially just observe a realization of a CM{ncx)N, p), whose typical clusters 
are GWp^p trees. 

All these examples show a stationary distribution for the typical clusters. However, 
this distribution is not stationary after gelation for any other example that we tried, but 
the computations are not as nice and do not make well-known distributions appear. 


A Index of notation 

Sets 

. N = {0,1,...}; 

. [iV] ={l,2,...,iV}; 

• M"'' = [0, -l-oo). 


Measures 

• Mi: probability measures on N; 

• M^^^: probability measures ^ on N such that Np{k) £ N for all A: £ N; 

• for p £ Ml, 

Gfiiy) = ^p{r)y'’, ?/£[0,l], 

r>0 

is the generating function of p] 

• for i £ N, 

< = G'«(1) 

are the factorial moments of p] 

• Fp{x) = xG'p{x); 

• for // £ Ml with > 0, 

/ ^ \ 
qii{x,y) = —X -hi- p\y; 


• for // / 0, 


p{k) 


{k + l)p{k + 1) 


A: £ N. 


Topologies 

• cadlag processes with values in a Polish space 5, endowed with the Sko- 
rokhod topology (see [7, 19]); 

• for /3 > 1, the space of sequences in N'^ with finite || • jj^-norm, where 

11^11/3= (1 + fcf+ + 
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for /c G N, C^: space of k times continuously differentiable functions on [0,1] or 
[0,1]^, endowed with the norm 


II/II7 


k 

sup 

i=0a:6[0,l] 


Dj{x) 


Smoluchowski equation We write 5 = N x N* and for p = (a, m),p' = (a', m') G S, 

• p ■ p' = aa'; 

• p o p' = [a + a' — 2,m + m'); 

• p' ^ P 0 ! < a + 1 and m' < m — 1; 

• if y ^ Pi p\p' = (a + 2 — a', m — m'). 

For two nonnegative families {c{p),p G S), {f{p),p G S), 

= IIc(p)/(p). 

p£S 

With a slight abuse of notation, (c, a) = (c, /) for /(a, m) = a, and so on. 


Discrete model 

• Np. number of particles in solution at t] 

• m = Nt/N: concentration in solution at t; 

• Pt{k,r): number of particles in solution k activated arms and degree r; 

• pt{k,r) = Pt{k,r)/N: concentration of particles that are in solution and with k 
activated arms and degree r; 

• ptir) = J2k=o Pt{k,r)/Np. distribution of the degree of the particles in solution; 

• TTt{k) = distribution of the number of activated arms of the par¬ 

ticles in solution. 

• ijjt- generating function of pp, 

• Mv) = i^tiliV) = 

• B{t) = J2k,r>0 kPt{k,r): number of activated arms in solution; 

• rp. i-th gelation time. 


Other quantities 

• CM(n,Tr^Gy_ configuration model on n vertices with degree distribution 

• pt = 1 - e“*; 

• qt{x,y) = {ptx + 1 -pt)y; 

• PN = N/a{NfG + {a{N)/N)^+^-, 

• the set of all couples 

(n, /i) G [J {n} X 

n>l 

such that 

cN < n < N, m 3 > m, '^2 r^^'^p{r) < M; 

r>0 
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set of (n, //) e k1^2,m such that additionally 


m 

m 


M 

1 


> 1 + m. 
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